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ABSTRACT 



Virtual Compton Scattering at High Energy 
ZHANG CHEN 

In this dissertation we develop a theoretical framework in the context of per- 
turbative QuantumChromoDynamics (pQCD) for studying non-forward scattering 
processes. In particular, we investigate a non-forward unequal mass virtual Compton 
scattering amplitude by performing the general operator product expansion (OPE) 
and the formal renormalization group (RG) analysis. 

We discuss the general tensorial decomposition of the amplitude to obtain the 
invariant amplitudes in the non-forward kinematic region. We study the OPE to 
identify the relevant operators and their reduced matrix elements, as well as the cor- 
responding Wilson coefficients. We find that the OPE now should be done in double 
moments with new moment variables. There are in the expansion new sets of leading 
twist operators which have overall derivatives. They mix under renormalization in 
a well-defined way. We compute the evolution kernels from which the anomalous 
dimensions for these operators can be extracted. We also obtain explicitly the lowest 
order Wilson coefficients. 

In the high energy limit we find the explicit form of the dominantly contributing 
anomalous dimensions. We are then able to solve the resulting renormalization group 
equations (RGE) and give a prediction of the high energy behavior of the invariant 
amplitudes. We find that it is the same as is indicated by the conventional double 
leading logarithmic analysis. 
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Introduction 



It is generally believed that Quantum Chromodynamics (QCD) is the fundamental 
theory for the strong interaction and it describes the properties of nuclear and subnu- 
clear, or hadronic, matter from first principles. QCD is a quantum field theory that 
is invariant under a local non-Abelian gauge symmetry-the SU(3) color symmetry. 
The basic matter particles of QCD, the quarks (and antiquarks), interact with eight 
massless gluons, non-Abelian gauge bosons that also couple to themselves, to form 
confined bound states of hadrons. The SU(3) gauge symmetry uniquely determines 
the forms of interactions between quarks and gluons, and dictates that the theory is 
both color confined and asymptotically free. This results in the bifurcation of QCD 
into Lattice QCD at low energy and Perturbative QCD (pQCD) at high energy. For 
more detailed discussion on various aspects of Quantum Field Theory and QCD, see, 
for example, References fl], ||, |], U, H, H> 0> H> HI • 

At low energies, because the coupling of QCD is strong, the usual analytic method 
of perturbative expansion in quantum field theory breaks down. As a result many 
fundamental properties of QCD at low energies, such as quark confinement, the dy- 
namical chiral symmetry breaking and ultimately the low energy hadron mass spec- 
trum, remain unresolved, at least analytically. The lattice formulation of QCD ]10 
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provides us with a new way of studying the theory non-perturbatively. By defining 
the quark and the gluon fields on a 4 dimensional Euclidean space-time lattice, the 
theory is regulated by an ultra-violet cut-off, the inverse lattice spacing. Physical 
quantities of interest defined through the path integral of the quantum fields can now 
be studied numerically, using Monte Carlo methods. 

At high energies, on the other hand, perturbation theory is valid due to the 
smallness of the QCD coupling. Perturbative QCD has had great success in explaining 
and predicting results of high energy experiments, for example, from the ratio of cross 
sections of e + e _ — > to e + e~ — ► hadrons, to structure functions of the nucleon 

in deeply inelastic scattering (DIS) (e.g., see [[II] and the references therein). In 
particular, pQCD has been successfully applied to the study of the structure of the 
nucleon, one of the most important frontiers in strong interaction physics. Here 
by structure we mean not only the traditional nuclear structure like spin structure 
and form factors but also, more relevant to the work in this thesis, the distribution 
functions of quarks and gluons, which we collectively call partons, inside the nucleon. 

Physics processes usually involve both hard (high energy) and soft (low energy) 
physics. The key to a pQCD analysis is the proper separation ("factorization") be- 
tween the hard part, where analytic calculation using perturbation theory is possible, 
and the soft part, where we have to rely on non-perturbative methods like numerical 
study and phenomenological modeling. By application of factorization theorems and 
resummation methods, and/or more formal methods like the operator product ex- 
pansion (OPE) and renormalization group equation (RGE) [|], ||, [12|, [Hfll , pQCD 
has been very successful in extracting the parton distribution functions of the nucleon 
and predicting their behaviors in a wide range of kinematic regions of high energy 
particle-particle scattering. However, there is still much work to be done, especially in 
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extreme kinematic regions like very small and/or very large Bjorken-x (see equation 
|1.7| ), where we have to not only look for new physics but also take greater care in 
separating the hard scattering process where pQCD is applicable and the soft physics 
parts where perturbation theory breaks down. 

Recently there is much interest in the analysis of non-forward scattering processes, 
for example, deeply virtual Compton scattering (DVCS) |15| , |l(| and hard diffractive 
electroproduction of vector mesons WA [13, [T^, ^| in DIS. This creates a new terri- 
tory to explore the quark and gluon structure of the nucleon besides the traditional 
inclusive (forward parton distributions) and exclusive (form factors) processes. A 
new collection of parton distribution functions called the skewed parton distributions 
(SPDs) [] appear in these exclusive, hard diffractive processes which contain new in- 
formation on long distance physics due to the non-forwardness (see, for example, Ref. 



2% P2|, |23| , |24|). At the same time, diffractive vector meson production also provides 



a new way of directly measuring the (forward) gluon density in the proton |l7j . The 
work of this thesis is aimed at developing a theoretical frame work to properly study 
non-forward processes in the context of operator production expansion and renormal- 
ization group analysis |25 |. The thrust of our research is obtaining the high energy 
behavior of these non-forward processes in general. 



DVCS and Skewed Parton Distributions 

The Compton process, referring to the elastic scattering of a photon off a charged 

object, historically provided one of the early evidences of the quantization and particle 

nature of the electromagnetic wave |2B| . Its role in studying the structure of hadrons 
1 SPD is the unified terminology used to replace the many terms like non-diagonal, non- forward, 
off-diagonal and off-forward from different parameterizations. 
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has been explored since the 50s with the derivation of Low's low-energy theorems 
27J which assert, for instance, that at sufficiently low energy the spin-dependent 



part of the Compton amplitude is determined by the anomalous magnetic moment 
of a composite system. Going to higher-order terms in the low-energy expansion, 
one finds the electric and magnetic polarizabilities. In recent years experimental and 
theoretical works in measuring and understanding the polarizabilities of the nucleon 



and pion have flourished f28fl . 

The virtual Compton amplitude for scattering a virtual photon off a hadron has 
become one of the basic tools in QCD to understand the short-distance behavior of 
the theory. The process of deeply virtual Compton scattering (DVCS) is, assuming 
the virtual photon is generated by inelastic lepton scattering, the Bjorken limit, i.e., 
the energy and momentum of the virtual photon going to infinity at the same rate, 
of virtual Compton scattering. The basic mechanism for DVCS is a quark absorb- 
ing the virtual photon, immediately radiating a real photon and falling back to the 
nucleon ground state. Ji |15|, ^] was the first to introduce and study DVCS. He pro- 
posed that we can obtain from DVCS information on the so-called off-forward parton 
distributions (OFPD) which in this case contain new information on long distance 
physics. 

The OFPDs have come up in different theoretical studies prior to Ji's work. In the 



late 1980s, Geyer and collaborators [29|, |3(J studied the relation between the Altarelli- 



Parisi evolution for parton distributions ]T3[ and the Brodsky-Lepage evolution |3T 



for leading-twist meson wave functions. The "interpolating functions" introduced in 



Ref. |29| are essentially Ji's OFPDs. In the early 1990s, Jain and Ralston studied 
hard processes involving hadron helicity flip in terms of an "off-diagonal transition 
amplitude" that involves off-forward matrix elements of bi-quark fields in the nucleon. 
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It was shown that the integral of this amplitude over the quark four-momentum 
yielded elastic form factors. 



Ji |21| introduced OFPDs in the study of the spin structure of the nucleon. 
His main observations were that the fractions of the spin carried by quarks and gluons 
can be determined from form factors of the QCD energy-momentum tensor, and that 
the latter can be extracted from the OFPDs. Furthermore, the DVCS process was 



proposed |15j as a practical way to measure the new distributions. Because the spin-2 
twist-two quark and gluon operators are part of the QCD energy-momentum tensor 
and because the form factors of the energy-momentum tensor contain information 
about the quark and gluon contributions to the nucleon spin, DVCS provides a novel 
way to measure the fraction of the nucleon spin carried by the quark orbital angular 
momentum, a subject of great current interest ||33|| . The more relevant part of his 
work to ours is his study on the evolution and sum rules (in terms of form factors) 
of OFPDs. He also obtained certain estimates at low energy of the OFPDs and form 
factors |22| using the MIT bag model [34fl . 



Radyushkin also studied the scaling limit of DVCS |16 |, and generalized the dis- 
cussion to hard exclusive electroproduction processes fl8], ^3[. The non-perturbative 
information is incorporated in his double distributions F(x, y; t) and non-forward dis- 
tribution functions F^(X;t). He discussed their spectral properties, the evolution 
equations they satisfy, their basic uses and some general aspects of factorization for 
hard exclusive processes. 

A third parametrization for the non-perturbative information was proposed by 
Collins et al |19], pCfl . 111 addition to the discussion of factorization theorem and evo- 
lution, they performed a numerical study of their non-diagonal parton distributions 
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in leading logarithmic approximation. They found that the non-diagonal gluon dis- 
tribution x%G(xi, X2, Q 2 ) can be well approximated at small x by the conventional 
(forward) gluon density xG(x,Q 2 ). 

It is clearly both interesting and important for us to study in depth these SPDs, 
because they will give us information on novel long distance physics due to the non- 
forwardness. 



Diffractive Vector Meson Production 

Conventionally, for example, in DIS, the gluon density xG(x, Q 2 ) is measured indi- 
rectly. The structure function F 2 (the same as Z/W2 of, eg, equation |1.91| ) that comes 
out directly from the experimental data is (essentially) the sum of quark (and anti- 
quark) distribution functions while the gluonic content can only be obtained from a 
DGLAP type evolution analysis (see Chap. 1 and Ref. jIT|). This is because gluons 
only have color charge and all probes in nature are color neutral due to color confine- 
ment. As a result any interaction with the gluonic content inside the nucleon has to 
have at least a quark loop involved and thus has to be of higher order. 

On the other hand, at very small Bjorken-x, the gluonic sector dominates the 
high energy behavior of scattering cross sections (see Chap.l for explanation). It is 
well-known that there is a sharp rise in the parton densities inside a nucleon at small 
x, which is driven by the gluons. While this might be an indication of the emergence 
of the so-called BFKL pomeron |35| , higher order DGLAP evolution is able to account 



for most of the increase (see, for example, Ref. [1 IJ). It is obvious that new ways of 
measuring the gluonic content in the nucleon, especially direct measurements, will be 
extremely useful. 
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It has been shown that the cross section for diffractive electroproduction of vector 
mesons can be predicted in pQCD [IT], |36], |37j (see Figure [I]). This process provides 



a novel probe of the dynamics of diffractive scattering in QCD. The prediction that 
the cross section is proportional to the square of the gluon density xG(x, Q 2 ) in the 
hadron [lTj means that we do have a process through which direct measurement of 



the gluon density is possible. Experimental data [[38] appears to be in accord with 
the predictions, including an enhancement due to the rapid rise of the gluon density 
at small x. 



The pQCD calculation leading to the above prediction in Ref. [|17j] was done under 
a double leading logarithmic approximation (DLLA)-resumming leading logarithms 
in both the longitudinal momentum (leading In 1/x) and the transverse momentum 
(leading lnQ 2 /A 2 , where A is the QCD scale). For further studies we need to relax (at 
least one of) these two approximations and also incorporate consistently the BFKL 
contribution. One such possibility can be found in applying the high energy factor- 
ization developed by Catani et al J}{|, relaxing the leading logarithmic approximation 
(LLA) on the transverse momentum k and discussing the cross section in terms of 
the so-called un-integrated gluon distribution (see, for example, [|0|]). However, up 



to the time of the work in this thesis |25[ , there had not been a completely consistent 
treatment. 

One of the basic problems lies in the fact that conventional forward factorization 
has been used on a non-forward process-because of the time-like four-momentum 
of the final state vector meson, diffractive vector meson production is necessarily 
a non-forward scattering process. It is therefore not obvious without theoretical 
investigation of non- forward factorization and/or OPE that it should be the forward 
gluon density that comes into the expression of the cross section. 
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The research of this thesis was indeed originally motivated by the diffractive 
vector meson production. The ultimate goal is a better understanding of the gluon 
distribution inside the nucleon, especially at small x. However, we realised that 
some formal theoretical development on OPE and renormalization group analysis in 
the non-forward cases has to be made before we can consistently attempt to relax 
the LLAs. In other words, we have to justify (or prove wrong) the use of forward 
factorization under DLLA before we can go further. 



Unequal Mass Virtual Compton Scattering 

This dissertation aims exactly at tackling the above problem. In this research, we 
formulate a general operator product expansion (OPE) description for a generic non- 
forward unequal mass virtual Compton scattering amplitude(see Figure 0). In order 
to avoid complications from dealing with meson light-cone wave functions (these have 
been shown by Collins et al |L9|| to be factorizable) , we have replaced the vector meson 



vertex with another virtual photon coupling. At this point, the main goal is to study 
the influence of the non-forwardness on OPE, rather than the phenomenology of the 
vector meson production itself. DVCS and DIS are extreme kinematic cases of this 
generic unequal mass double virtual Compton scattering (see section |2.1|) . 

General proofs of the DVCS factorization have been given by Radyushkin in his 



approach based on a-representation [16, 18, E3l and by Ji in an alternative method 



nj. It has been established that factorization for DVCS is on the same footing as 
other well-known examples like DIS. Collins et al [[HJ have also proved the factoriza- 
tion theorem for general diffractive vector meson production. 

Some early studies of unequal mass Compton processes can be found in Refs. 
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|. The factorization property of the general two virtual photon process can 
be summarized beautifully in terms of Wilson's OPE. This expansion now requires 



operators with total derivatives |25], |3(| |4!| to describe the non-forward nature of 
the process. While it is well-known that these derivative operators contribute to the 



wave functions of mesons |3T], |43|, and that the scale evolution of total derivative 
operators can best be studied using conformally-symmetric operators our work 
will concentrate on the renormalization group evolutions of these operators (and their 
corresponding Wilson coefficients) in the high energy (and small x) limit (see later 
chapters of this thesis). 

Because of the non zero momentum transfer r = p — p' from the initial state 



proton to the final state proton (Fig. [Tj^), the scaling behavior must be different 
from that in a forward case because of the new kinematic degree of freedom (the non- 
forwardness) . Indeed we will show that compared with forward scattering one has two 
scaling variables, or equivalently, two moment variables u and v (see equation 3.13| ). 
The amplitude should now be expanded in terms of double moments with respect to 
these two moment variables. There are associated with these double moments new 



sets of operators that have overall derivatives in front (see equation |3.11|) , which mix 
among themselves under renormalization via an anomalous dimension matrix. The 
reduced non-forward matrix elements of these operators are thus double distribution 
functions. However, they do not seem to have a simple probability interpretation. 
While formally we can go from DVCS to DIS using the same formalism and in principle 
we can proceed beyond DLL to next to leading order, we will focus our attention on the 
high energy behavior of these distribution functions. We solve their renormalization 
group equations at high energy, and show that in the high energy limit the gluonic 
double distribution function actually reduces to the conventional (forward) gluon 
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density and the high energy behavior is the same as obtained by a conventional 
double leading logarithmic (DLL) analysis. This also means that we have justified 
the DLLA calculations JHj| . 

The outline of the thesis is as follows. In Chapter 1 we review the standard 
operator product expansion and renormalization group equation analysis of forward 
scattering amplitudes, using DIS as an example. In Chapter 2 we define the process 
and the kinematics of the non- forward scattering amplitude T Mi ,. A general tensorial 
decomposition of T^ into invariant amplitudes is given. In Chapter 3 we perform 
a general operator product expansion of T MI/ , define new moment variables and give 
their relationship to the more conventional Bjorken type scaling variables. In Chap- 
ter 4 we write down the renormalization group equations, present the equivalence of 
evolution kernels for the double moments, and calculate explicitly the lowest order 
Wilson Coefficients. In Chapter 5 we go to the high energy limit to solve the renor- 
malization group equations and make connection with the double leading logarithmic 
analysis. Chapter 6 gives conclusions and outlook for future work. 
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Chapter 1 
Forward Scattering 

In this chapter, we give a review of the traditional operator product expansion (OPE) 
and renormalization group (RG) analysis for a deeply inelastic scattering (DIS) of 
lepton and proton (e-p), which serves to lay out the framework and conventions that 
will be used in later discussions of non-forward scatterings. This will be the forward 
case that we refer and compare to in those chapters. 

We start with the parton model description of DIS to present the problem and 
establish the kinematics. We proceed to decompose the amplitude into invariant 
amplitudes and define the structure functions. We then perform an OPE of the 
invariant amplitude and discuss the RG equation and its solution of the operators 
and Wilson coefficients involved. To set the conventions, we also include a review 
of the Light Cone(LC) gauge convention and explicitly calculate the quark-quark 
splitting function in the LC gauge. 

Most results in this chapter concerning quantum field theory are more or less standard 
material and can be found in, for example, 0, [|j] or ||, while more specific details of 
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DIS process can be found in for example, |ll|. However, we will be presenting them 
in a different way ||45|| . We will generally not make explicit references to the standard 
literature on well-known results. 



1.1 Amplitude and Kinematics 

The process we consider here is an electron-proton scattering at very high energy, 
such as those that occur at DESY. Figure || depicts one such process, where the 
initial proton has four momentum p and spin a while that of the initial electron is 
k and A, respectively. After the scattering, the electron has four momentum k' and 
spin A' while the final state of the proton is collectively labeled n. 

The exchanged virtual particle, with momentum q = k — k', can be either a photon, 
a Z or a W boson (in figure ^ it is a photon). For the purpose of our discussion, we 
will consider only one photon exchange, which is a good approximation at the leading 
level. In the kinematic region we will be working, that is, with Q 2 = —q 2 > 1 GeV 2 
and energy of the electron at tens of GeV (for example, DESY has a electron beam 
energy of about 30 GeV), we can neglect the mass of the electron in the following 
discussion, which means k 2 = k' 2 = 0. 



1.1.1 Invariants of DIS 

We first work in the rest frame of the proton to define the kinematic variables used 
to describe such a process. We have 

U=(M, 0,0,0) 
\ k = (e, 0,0, e) 



13 



If the scattering angle of the electron is 9, then 

fc = ( e', e' sin 6*, 0, e' cos 6*) , (1.2) 

thus the exchanged four momentum square is 

Q 2 = -q 2 = —(k - k') 2 = 2k ■ k' - k 2 - k! 2 = 2ee'(l - cos 9) . (1.3) 

The first invariant of a DIS process, Q 2 , is then given by 

Q 2 = 4ee' sin 2 °- , (1.4) 

which indicates the hardness of the scattering. We can see that large Q 2 needs large 
e, e', and 6 1 not approaching zero. 

A second invariant v is defined by 

Mu = p-q = P -(k-k') = M(e-e), (1.5) 

which gives 

z/ = e-e'. (1.6) 

We can see that v is the energy loss of the electron in the proton rest frame. 

A third invariant, the Bjorken-x variable, is more commonly used instead of v 



and is defined as 



We note < x B j < 1- 



= 1 Q^ = Jl_ h7 , 

XBj ~ u~ 2Mu~ 2p-q- [ } 



Q 2 and XBj are the usual variables used to describe DIS, and sometimes we use 



equivalently Q and v. Although they are defined above in the photon rest frame, 
they are by definition invariant when we go to any other frame. 
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1.1.2 Scattering Amplitude and Differential Cross-section 

The scattering amplitude of the DIS process in figure can be written as 

T^ A = ^U(^U{k)^^{n^\j v {0)\pa)ie, (1.8) 

with the minus sign on the n state meaning outgoing states and j v the electromagnetic 
current given by 

3» = J2 e fQtf»Qf- (1-9) 
/ 

The differential cross-section is then given by 



d° = ~ EE |7SaI 2 (2-) 4 ^ + q - Vn)^ , (1.10) 



A' A o~n z,tz,t 



where the factor | comes from averaging over initial proton and electron spin and the 
phase space and kinematic factor are obtained from a plane wave normalization. 



To simplify the expression of da we define the electron factor and the proton 
factor such that 

dk« ~ 2e2e> q* L ^ ■ I 1 - 11 ) 

We have, for the electron part, 

U = ^E^(^)7^a(A:)]*[^(^)7^a(A:)] 

/ A'A 

= lj2Ux(kh fl U x/ (k , )U y (k') lu Ux(k) . (1.12) 

^ A'A 

Because we know that 

E U X ,a(k)UxAk) = {¥ + ™e)a/3 , (1.13) 
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we arrive at 



1 / 
lu.u = 2 Tr ^(^' + m e )>y u (f{ + m e )} 

= Irr[ 7 ^'7, (1-14) 



where we have taken the limit m e = 0. Thus we obtain 



lf„, = 2(k'^k u + k^k' v - g^k ■ k') + 0{m 2 e ) , (1.15) 

which is determined completely from the kinematics of the process. 
On the other hand, for the proton part, we have 

ir = ^EE{^ ( " ) b' /i ( )b (7 )*^ ( " ) l/( )l^)( 27r ) 454 (p + ?-Pn)} 

^ a n 

= ^E^I^WI^^^I^WMls^b + ^-Pn), (i.i6) 



where we have used the fact that j M is Hermitian. 

Because the presence of the 5-function, the sum of n state is not a complete set 
summation, but rather an on-shell summation. We rewrite the (^-function and obtain 

IT = Wf d 4 xe^ +q -^>^ P a\f(0)\n^)(n^\f(0)\pa) 

= III I d 4 xe^(pa\f(x)\n^)(n^\f(0)\pa), (1.17) 

where we have used the translation relationship 

j,(x) = e *%(0)e-** . (1.18) 

This translation operation is exactly analogous to the spacial translation by the three 
momentum operator p or the time translation by the Hamiltonian H in quantum 
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mechanics (see, for example, |[46| ) . In the above, p is the 4-momentum operator 



where 



p ■ x = Ht — p ■ x 



\ p\ Vacuum ) = 
^ Pti\n) = Pnn\n) ■ 



(1.19) 



The n state summation is now of a complete set, thus by the completeness relationship 
we arrive at the expression of the proton factor as 



1.2 Structure Function 

The structure functions of the proton are defined from the invariant amplitudes of 
the proton factor, which in turn are obtained by a general tensorial decomposition. 

1.2.1 Tensorial Decomposition 

L^v as a rank-2 tensor depends on the kinematic variables p and q only. Thus the 
general tensorial decomposition of L^ u can be written as 



L^u(p, q) = + Bp^p v + C{p^q u + p u q^) + Dq^q u + E{p ti q v - p v q^) + Fe^ vpa q p p a , 

(1.21) 



where A,B,... are invariant amplitudes depending only on the invariants of the pro- 
cess, that is, A = A(p ■ q, q 2 ) and so on. We have deliberately grouped the terms into 
symmetric and anti-symmetric parts to make the symmetry properties of L[iv more 
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manifest. The last term does not occur for an electron scattering, however, we do 
have to consider it if it is a neutrino scattering. 



The conservation of the electromagnetic current gives (cf equation |2.26| ) 



<f L„„ = q v L™ = . 



1.22) 



Thus we have 

<f V = Aq v + Bp ■ qp v + C(p ■ qq u + q 2 p v ) + Dq 2 q v + E(p ■ qq v - q 2 p u ) = , (1.23) 

being true for any p and q. We can vary p and q while keeping p ■ q and q 2 constant 
and this means that the coefficients of the individual p v and q v have to vanish. We 
have 

A + Cp ■ q + Ep ■ q + Dq 2 = 
Bp-q + Cq 2 - Eq 2 = 

Similarly, from 



1.24) 



q v L ilv = Aq^ + Bp ■ qp^ + C(p ■ qq^ + q 2 p il ) + Dq 2 q ll + E(q 2 p^ - p ■ qq^) = (1.25) 



we have 



A + Cp ■ q - Ep ■ q + Dq 2 = 
Bp-q + Cq 2 + Eq 2 = 



1.26) 



Combine equations |1.24j and |1.26| we arrive at 



E = 

A + Cp ■ q + Dq 2 = 
Bp- q + Cq 2 = 



1.27) 
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which eventually leads to 



C = - 
D = - 
E = 



B) 



:i.28) 



Therefore, only two of the invariant amplitudes survive after considering current 
conservation and we have 

P ' Q 1 (p * <?) 2 
Lfj,u = Ag^ u + Bp^p v - B—^-{p il q v + p u q^) -{A — — B)q fl q v , (1.29) 



or eventually, after collecting terms, 



V = Afe. - q -f) + B - + m.) + fef ft* 



q A 



1.30) 



1.2.2 The Structure Functions 



The traditional structure structure function of the proton is indeed just L^ v with 



kinematic factors. Explicitly (see 1.20Q 
1 M 



W^(p,q) =L, v = l -Y,j d 4 xe^(pa\nx)f(0)\pa) . (1.31) 



And after the tensorial decomposition, the two proton structure functions Wi and 
W2, which are only functions of the kinematic invariants of the problem, are defined 
by 



W 2 (Q 2 ,x) ( 
+ M 2 W 



q 2 

P-Q 

Q 2 



(9 



1.32) 
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Now in terms of the structure functions, the differential cross-section |1 . 1 1| becomes 

d 3 a 1 M 1 e 4 



W Mi 7 

dk' 3 4tt 2 E p 2e2e' Q 4 ^ ' 



1.33) 



In calculating W^l^ we note that q^ v = + 0(m 2 ), so we can drop all terms 
having explicit factor of or q u in W^. We get 

Wo 



4WiA; ■ k' + 2 jfi( 2 P -kp-k' - M 2 k ■ k') 

6 6 
4ee / (2Wi sin 2 - + W 2 cos 2 -) . 



1.34) 



Thus the differential cross-section is given by 

d 3 a 



A o> 2 M ,„ TTT .08 „ T 9 6, 
—— = — (2Wi sin 2 - + W 2 cos 2 - 
dk' 3 Q 4 E P K 2 2' 



1.35) 



where a em = |- is the fine structure constant. 



It is conventional to give da jdQ 2 dv rather than da/d 3 k' , so we need the Jacobian 
of the variable transformation. We have 



d 3 k' 



27rdcos9(e'yde' 

d cos 9 d cos 9 



2n(e 



l\2 



dQ 2 

cV 
dQ 2 



du 

d£ 
du 



dvdQ 1 = Ti^ydvdQ 



dQ 2 


dQ 2 


d cos 9 


cV 


du 


du 


d cos 9 


cV 



-1 



Tr—dudQ 2 . 
e 



1.36) 



Therefore, the expression of the differential cross-section eventually becomes 
da a a lm e ' M ,„ut • 2 9 ttt 6\ 

W = 4l #W 2W > sm V W2COS V 



1.37) 
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In the laboratory frame where 6 is the scattering angle, it is possible to fix x and 
Q 2 (or equivalently Q 2 and v) and vary 9, thus experimentally Wi and W2 can be 
determined separately. Since there are more small angle experimental data, W2 is 
better measured. 

Theoretically, what we usually do is to relate the cross-section, in terms of the 
structure functions, to the imaginary part of a forward scattering amplitude |L2] via 
an optical theorem. 

To be more specific, define the forward scattering amplitude T^ u as 

T, v = *4vr 2 |^ £ / d 4 xe^(pa\Tf(x)nO)\P<?) , (1-38) 



we have (see |1.31| ) 



W p = 2/mT F . (1.39) 

The forward amplitude is shown in figure [3](a) while the structure function is shown in 
figure 0(b), both up to kinematic factors. The only difference in the graphs is the cut 
in the middle in figure 0(b), which means that all the intermediate lines the cut goes 
through are put on the mass shell. Mathematically this is equivalent to replacing the 
Feynman propagators of all the intermediate particles with its corresponding on-shell 
^-function and the proper factors (see section |L4.3|) . 



1.3 Operator Product Expansion Analysis 

We will use a straight forward operator product expansion (OPE) analysis to discuss 
the deeply inelastic scattering (DIS) of previous sections. While similar and more 
complete discussions can be found in H, we will present ours in a slightly different 



notation, adopted from [45 . 
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1.3.1 General Statement 

Let 0\(x) and 02(x) be local operators built out of fundamental fields like q, q (quark 
and anti-quark fields) or (gauge fields) and finite number of (covariant) derivatives 
at one space-time point, for example, 

jfa) = q f {x)l^q{x) ; d^(x) = qD^(x)q{x) ; ... . (1.40) 

We have, in the short distance limit, 

6 1 (x)d 2 (0) ^EC.,J ) Ejr^Jx) + Remainder . (1.41) 

r=l 

are local operators which can include the identity operator I and E^ r \x) are 
the so-called Wilson coefficients, c-number functions that are usually singular in the 
limit of x — > 0. 

By a dimensional analysis we can obtain the asymptotic behavior of the Wilson 
coefficients E^'\x) at short distance. The convention we use is such that mass and 
energy has dimension +1, i.e., [E] = dim(E) = [M] = dim(M) = +1, while co- 
ordinate has dimension of — 1, ie, [x] — — 1. Thus, we have under our convention, 
[q] = — |, [Ap] = —1, and so on. 

Now let d r be the naive dimension of 0^ r \ di and d 2 be those of Oi and 2 , respec- 
tively, E^ r \x) has an x behavior/dependence, for small x, given by (v / a?) _A , where 
A = d r — di — d 2 since we must have d\ + d 2 = d r — A. That is, we have 

E tL-,J x ) X - (Vx~ 2 ) dl+d2 ' dr ■ (1-42) 

When we increase N, we can only generate new operators by adding new fields or 
more derivatives into the operators, these new operators must have more negative 
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dimensions, that is, d r decreases with increasing N, therefore, we may choose N large 
enough such that the remainder becomes as small as desired, namely, remainder ~ 
(x 2 ) R for any desired positive value R in the small x limit. 



1.3.2 Renormalization Group Equations 

We work in the zero quark mass limit. The operators under renormalization behaves 
like 

Oi(x, li 2 , = Z i \— , a M0 ) O^x, /xg, a Mo ) , (1.43) 
Mo 

where we have explicitly shown the dependence of the operators on the renormaliza- 
tion scale. This is the so-called multiplicative renormalization and Zi more generally 
can be a matrix when there is mixing among the operators. 



Taking the logarithmic derivative of the scale fi on both sides of 1.43 we have 



1 2 dZi -y * . 2 \ 

li(<*i) di(x,fi 2 ,ap) , (1.44) 



where in the last step we have defined the so-called anomalous dimension 7^ of the 
operator O, as 

■*m = (1 - 45) 



Note that 7, is a function of a M only. This is because from |1.45| 7« can only be a 
function of fi 2 and ct M and is dimensionless. Because we have taken all the quark 
masses to be zero, there is no scale to set (/,, thus 7^ is a function of the coupling only. 

2 

On the other hand, if m q 7^ 0, then 7^ — > 7i(a M , ^-2), complications will arise. 
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The operators, at the same time, obey a renormalization group equation of the 
form 

At 2 ^Oi(x, fx 2 , a M ) = 7i(a Ai ) Oj(a;, ,u 2 , a„) , (1.46) 
or more generally, when there is operator mixing, 

V 2 -^di(x,fx 2 ,a^) = Y^7ij(®n) dj(x,fx 2 ,a^ . (1.47) 

To establish the renormalization group equations of the Wilson coefficients, let us 
look at a simple case where only one term is kept in the operator product expansion, 
that is, 

6i(x,/x 2 ,a M )d 2 (0,/x 2 ,Q; M ) 2=3 6 {3) (0,/i 2 ,a M ) E (3) {x, fi 2 , a J + Remainder . (1.48) 

Taking the logarithmic derivative of /z 2 on both sides and omitting the contribution 
from the remainder, we have, suppressing the /i 2 and dependence, 

(7i + l2)d 1 (x)O 2 (0) ^ l3 6^(0)E^ + 6(3) (0)^-1^(3) . (1 . 49) 



Combined with the operator product expansion itself, we have 

^ 2 ^E^(x) = {^A- 2+ ^)E^{x) = (7i + 72 - 73)^ (3) (^) , (1-50) 
where (3 = P(a^) = /J 2 ^ is the QCD beta-function. 

Therefore, the renormalization group equation for the Wilson coefficient functions is 

(d d \ 

/^ 2 ^ + /5^ + 73K) -7iK) -72K) ) EV)(x,n 2 , aiM ) = 0. (1.51) 
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1.3.3 OPE Description of DIS 

We now put the operator product expansion formalism to use on the deeply inelastic 
scattering process in Q| We work on the forward amplitude defined in 1.38 and 



the structure functions (and thus the cross-sections) can be obtained through |1.39 
As usual we work in the so-called Bjorken limit where we have 

— q 2 = Q 2 large 



mv = p ■ q 



large (1.52) 



%Bi = — — o fixed . 
Forward Operator Product Expansion 

We expand Tj Al (a;)j !y (0) as x^ — > 0. The tensor structure of T Mi/ can be proven as (|5 
Tj»(x)j u (0) = A, p E^\x 2 )I + F%\x 2 )6^}(0) 



) *Xj lAj a a a *Xj 



+ B^Y: E^\x 2 )O% n l^ n (0)x^x^...x^ + Rem. , (1.53) 

i,n 

where A^ u and are conserved tensor structure operators and (see also section 



B^vap = g m Qvfi u + g^udad/3 - g m d v d p - g v pd^d a . (1.54) 

We suppose symmetric combinations in /i, v (since from section |1.2.1| , T^, after all, 
is symmetric) and also we suppose that all indices in Os are symmetrized (even a, 
(3 with /ijs). We note that the Os are the same operator sets for both terms and 
that the label n is actually the angular momentum quantum number, or spin, of the 
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corresponding operator. The term twist is defined as the difference between the naive 
dimension of O and its spin, ie, twist of O n = [O n ] — n. 

We need to evaluate the above OPE between two symmetric external proton 
states (see |1.38|) , however, in practice in the kinematic region where high energy 



particle experiments are conducted, especially in the small xbj region where most 
of our interest lies, we can put in some transverse momentum to make the external 
states asymmetric and physics will not see the difference at this level (note it is the 
main goal of this thesis to discuss what will happen with asymmetric external states 
when we go to higher level of accuracy). That is, 

{p\Tjp(x)j„(0)\p) = \im(p - r|Tj M (z)^(0)|p) . (1.55) 
Therefore the identity operator I does not contribute to in the OPE. We need to 



work to the next terms in 1.53 and we have 



i,n 

+ (g^a9^ a + g^udad/3 - g m d v d p - gvpd^da) 

£ E£\x 2 )x^x»\..x^( P \dt^Sm ■ (1-56) 

i,n 

For the matrix elements {p\O^ n ^ n (0)\p) , the indices can be made only from either 
g^ u or because x M is gone, however, the contribution of the two types of indices are 
different. After the Fourier transformation, we know from dimensional analysis that 



q 

Thus we have, for the terms, 



%. (1.57) 
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while at the same time, 

qz gz qz 

Therefore the g^ u contributions are small compared with those of p^p v . So the 
only term we keep is the contribution from a totally symmetric (in //j) combina- 
tion p^p^-'-PfMn- ^ is clear that this approximation is good to O(^). That is, we 
write 

(P|0£i„(0)b> = TVJW^ (P\\0^\\p) + 0(1) , (1.60) 

where (p||OW||p) is the so-called reduced matrix element of the operator O^^. 
Because the external states are on-shell protons, (p||Ow||p) does not depend on kine- 
matic variables and thus is a number depending only on i and n, the flavor and spin 
indices, respectively 

Note indeed this approximation is the so-called leading twist approximation be- 
cause non-leading twist operators will have in their contributions extra factors of ^ 
following similar discussion as the above. 

We therefore obtain, after taking leading twist approximation, 

(p\Tj,(x) Ju (0)\p) ^ £ {{g, v U - d,d u ){p-xTF^{x 2 ) 

i,n 

+ [g,u(p-d) 2 +p^ Pu a - (p^d u +p v d^)p-d](p-x) n - 2 Etl 2 (x 2 )} . (1.61) 

The forward amplitude T^ is essentially the Fourier transform with momentum 
q of the matrix elements (|1.38|) . And under the Fourier transform, it is clear that 
x^ =5- Thus we can rewrite the part concerning the Fourier transform of the 

Wilson coefficients E^(x 2 ) in terms of logarithmic derivatives of q 2 as 

[ d*xe iq - x (p-x) n E%\x 2 ) = (-ip^Y f d 4 xe iq - x E^(x 2 ) 
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(H2p-g)A r g«(^ 



-i2p ■ q\ ( 2 9 



n 



\ ,r J \ q W) 4 (L62) 

where is the Fourier transform of E$\ and we have a similar set of equations for 
FW. (Note in our notation, (g 2 ^)™e«(g 2 ) is exactly the coefficient (7 Jo (Q) in § 
with the replacement of n — > J (spin index) and % — > a(flavor index). ) 

We then define the Wilson coefficients in the momentum space E$ and as 

J d^xe^ip-xTF^) = ± (^YpWffl ee ± U »F®ffl , (1.63) 



where we have used equation ET , Substitute these definitions into the Fourier trans- 
form of the matrix elements and note that using integration by parts we can show that 
the derivatives in the conserved tensor operators ( |1.54 ) become factors of momentum 



q in the fashion <9 M — > iq^, we arrive at 

/ d'xe^( P \TUx)Uo)\p) = -Y,i(p\\d®\\p) {- (g»v - ^r) ^f^\q 2 ) 

+ I ~ — ^ ~P ' 9 + 9^—^— J w E K n l 2 (Q ) | . (1.64) 

Rewriting the last term in the parenthesis in front of E^} 2 as 

»^ = (»,.- + (1-65) 

we obtain 

d'xe">*{ I j\Tj l Xx)jM\p) 



2\ 



+ \PnPv ~ 2 P ■ Q + — 2 ^- 1 ^ ^(Q ) j • (1-66) 
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From equation |1.7| we have 



2 (p-q) 2 1 



Q 2 u 2 q 2 2 ' 



1.67) 



thus 



d 4 xe^(p\Tj,(x)j u (0)\p) 

( K I'A L, 



-f\ u n {F«\Q 2 )+\E%m) 

tt f ^ - ^(p^+P^) + T^q,qu) u^ l E%{Q 2 )) (1.68) 



Mz/ 



(5 



Recall the definition of T„^, equation 1.38, and the tensorial decomposition of 



the structure function, equation |1.32| , we obtain the tensorial decomposition of the 
forward amplitude as 



1.69) 



where the forward invariant amplitude Ti^ are related to the structure functions Wi,2 
by 

Wi = 2/mTi ; W 2 = 2/mT 2 . (1.70) 



Comparing with the expression of the operator product expansion of T^ u , equation 
|1.68| , we obtain the formula for the invariant amplitudes as 



Ti 



uT 2 = ^E^^T.iPWO^Wp) B&(Q 



1.71) 
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Dispersion Relation and Optical Theorem 

Equation |I.71| is essentially a power series expansion of the invariant amplitudes. To 



further the computation we need the analyticity properties of Ti and T 2 

Explicitly write out the time ordered product in T^ (suppressing the spin aver- 
age) by inserting a summation over arbitrary intermediate states |r) we get 

V = tAn 2 ^Jd 4 xe^(p\Tf(x)f(0)\p) 

= 47r2 fp£ / d 4 xe^{ 6(x )(p\ 3fl (x)\r}(r\j u (0)\p) 

+ 9(-x )(p\j„(0)\r)(r\Mx)\p) } . (1.72) 



By applying the translation operators (see equation ETj ) we have 

V = ^vr 2 ^^ J rf 4 xe^{^(xo)e*-^-(p|j,(0)|r)(r|j,(0)b) 

+ e(-x )e-^-^(p\UO)\r)(r\j,(0)\p) } 
= ^7r 2 ^(2vr) 3 ^ J ^ {5 3 (g>p-p r )e l ^^^^(xo)(p|j M (0)|r)(r|j^O)b) 
+ 5 3 (q -^+p r )e i ^-^°)- a: ^(-xo){pb,(0)k}(rb M (0)b) } 

= i^ Ep (2^) 3 Y ( 6 ^+P-Pr)(p\3M\r)(r\]M\p) 
M , \ g + Po - Pr,o + ie 

_ S\q-p + Pr)(p\]umr)(r\]M\p) \ (1 ?3) 

go -po+Pr,o -*e / 

where in the second step we have performed the spacial integration that yields the 
3-d momentum 5-function. In the last step when we perform the time integration we 
have put in the proper damping factors ±e in the exponent to assure the convergence 
at the boundaries in an adiabatic approximation. 

To visualize the analytic properties of T^, again we go to the rest frame of the 
proton where p ■ q = Mq = Mv. The independent variables T^j, depends on are p 
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and q while for the invariant amplitudes the independent variables are actually Q 2 
and uj, with now (see equation |1.7|) 

2p-q 2M 

Thus, for fixed Q 2 , singularities in uj are the same as singularities in q . Note this 
result is covariant-it is just for the purpose of visualization that we need to choose a 
frame. 



From equation |1.73| we have singularities in q at (i) qo + po = p r ,o = E r and 



(ii) qo — Po = —E r . We shall discuss them one by one. 
(i) The first set of poles in q occur when 

(g + Po? = E 2 r =M 2 +f r = M 2 r + (q + p) 2 , (1.75) 

where in the last step we have used the spacial (^-function. M r is the invariant 
mass of the intermediate state |r). We therefore have 

M 2 r ={q + p) 2 = -Q 2 + 2p ■ q + M 2 , (1.76) 
and the corresponding poles in uj are 



Ml - M 2 

~Q 2 



u>r = l+ r n2 ■ (1-77) 



(ii) The second set of poles in q Q occur when, after similar discussion as above, 

M 2 = (q-p) 2 = q 2 -2p- q + M 2 , (1.78) 

and the uj poles are at 

Ml — M 2 , 

Ur=~±- r Q2 • (1-79) 
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The intermediate states |r) can be \p), \p,ir), \p,n,7i), etc. Electro-magnetic current 
and Baryon number conservation require that M, mass of the proton state), be 
the smallest invariant mass of all possible states involving any baryon. This means 
we always have M r > M. Therefore we can clearly visualize the analytic property of 
T Ml/ on the complex cu-plane as shown in figure |j: 

• has a pole at u — 1 and another one at u — — 1. They correspond to an 
elastic scattering where |r) = \p) and thus M r = M. 

• also has two branch cuts. One on the positive real cu-axis and starts from 
the point to = 1 + ( M+m ^) ~ M ; which corresponds to the lowest excited state 
\p, 7r) with a pion generated almost at rest and as such M r = M + m^. It 
extends out to infinity, or rather, as long as the collision center of mass energy 
is big enough to generate the states. The other one is simply a mirror reflection 
of this cut about the imaginary cu-axis. 

It is therefore obvious that T^ has an analytic circle of unit radius on the tu-plane 
and thus its Taylor expansion in uj exists as long as we are inside the unit circle. 

However, the region of < u < 1 is not physical. The physics, namely, deeply 
inelastic scattering in the Bjorken limit, is happening in the kinematic region where 
< %Bj < 1 (see section |1.1| ) and corresponds to uj > 1. Thus the operator product 
expansion of equation |1.71| , although called the OPE of the forward amplitude, can 



not be applied directly as it is to the physical scattering, nor the computation of 
structure functions and cross sections. 



It turns out |12j that by deriving a dispersion relation for Z/T2, we can indeed 
correctly apply the operator product expansion to the physical situation and relate 
directly to the structure function 
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We start by writing a dispersion integral of vT 2l 

uT 2 (u;,Q 2 ) = -L / -^zyT 2 (^,Q 2 ) (1.80) 
2m JcOJ — lo 

where c is any coutour enclosing the origin on cu-plane and lying complete inside the 
unit circle (see figure f|). The analyticity of T M „ within the unit circle will guarantee 
its convergence. 

We can continuously and analytically distort the contour c to d as shown in figure 
H|, where d crosses no cut, encloses no pole and has its boundary segments pushed 
to infinity. At infinity although actually z/T 2 approaches a constant [jj5|, there are 
cancellations between ±00 resulting in extra convergent factors and as such we can 
drop the integral region of the two semi-circles at infinity and obtain 

vT 2 (u,Q 2 ) = -^T 2 (u/,Q 2 ) 

Z7U Jc->c' U) — U) 
I roo A, / 

= — (uT 2 (Lu' + te,Q 2 )-uT 2 (iu'-te,Q 2 )) 

2m J 1- uo — uo 

1 r— OO A, / 

+ — / -=— (*T a (a/ ~ fe, Q 2 ) ~ ^ 2 {J + ie, Q 2 )) ,(1.81) 
2m J-i+ 00 — lo 

where e is a infinitesimal positive number. 

To simplify the above expression, we need the symmetry property of T Mi , in uj, or 
equivalently, in q . as a time-ordered product can be represented systematically 
by a set of Feynman diagrams. Equation |L73j is for real go values. If we take go to 
be complex and forget ie, we have an extended T^g^) on a complex plane. The 
usual time-ordered product is obviously obtained when we take go — > I go I + it or 
Qo — I go I _ it- This extended T^(gg) clearly obeys 

V(g ) = V(-go) , (1.82) 
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which can be obtained directly from equation 1.73 . When we go to real values, we 
actually have 

T 2 {uj + te,Q 2 ) = T 2 {-u-te,Q 2 ), (1.83) 



which means that for z/To we have 



uT 2 (uj ± ie, Q 2 ) = —uT 2 (—u =F ie, 



1.84) 



Therefore, by setting to' — > —uj' in the second integral of equation |1.81| and use the 
above we obtain 



vT 2 (uj, Q" 2 



1 f°° du' 



2m Ji- uj' — uj 

1 r°° du' 



(z/T 2 (u/ + ie, Q 2 ) - vT 2 {u' - ie, Q 2 )) 



ivT 2 (-u' - it, Q 2 ) - vT 2 (-uj' + ie, Q 2 )) 



+ n • / 

2m Ji- uj + uj 

1 roo A, / 

' {uT 2 {u' + ie,Q 2 )-vT 2 (uj' -ie,Q 2 )) 



-vT 2 (uj> + ie, Q 2 ) + vT 2 (uj' - ie, Q 2 )) 



2m Ji- uj' — uj 

1 r°° du' 



+ 



2m' Ji- uj' + uj 

JL r d J (__! _L_ \ [uT 2 ( U j'+ i e,Q 2 ) - uT 2 (uj'^e,Q 2 )} . 

2m J \- \u'—uj uj'+ujJ 

1 /-co 2lj 

— / dJ— -\v^ 2 (uj'+i^Q 2 ) - vT* 2 (uj' + te,Q 2 )] 

Z7TZ JI - Cl> — UJ 

uj r 00 duj' 



m J i- uj' 2 — uj 2 



2i Im vT 2 (uj' + i eps, Q 2 ) 



1.85) 



where in the last step we have used the Hermiticity of the electro-magnetic current. 
Recalling equation |1.39| , the optical theorem, we obtain the dispersion relation of uT 2 

as 

i,i roo A,,) 1 

(1.86) 



uT 2 (uj,Q 2 ) = - r -^- 2 uW 2 (uj' } Q 2 ) . 

IT J I- UJ Z — UJ Z 



Note that in the above, uj should still be seen as inside the convergence circle (the 
unit circle) of uT 2 while only uj' has the physical meaning of the inverse of Bjorken-x. 
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Therefore we can expand the denominator of the dispersion integrand and get 

uj r°° duo' „, , , „„ 1 



uT 2 (uj,Q 2 ) = - r%uW 2 (u',Q 2 )- 

7T J I- UJ 1 



uj [°° duj' Txr ^ /cu 2 \ n 



7TJ1- u/ 2 ^oW J 

E — / T^^KQ 2 ). (1-87) 



Change the integration variable from a/ to = u' = x we have 

do;' 



dx 



0J> 2 



iU x 2n dx, (11 



(uj') 2r * 2 



and thus 



uT 2 (uj,Q 2 ) "= (0) fv^i [ 1 dxx 2n pW 2 (x,Q 2 ) 

n=0 

oo i ! 

^ uj n+1 - dxx n uW 2 (x,Q 2 ). (1.89) 



n even 



Again we would like to stress that uj and x are completely different quantities, uj 
is essentially a mathematical qunatity introduced to Taylor expand the (invariant) 
amplitude(s) around the origin in an operator product expansion. It is the so-called 
moment variable and has its value limited between and 1. On the other hand, x 



is the Bjorken kinematic variable of the actual physical scattering (see |1.7| ) 

with its value also limited between and 1. The convolution of the structure function 
uW 2 with the n + l-th power of x is called taking the n-th moment of the structure 
function. We can see that because of the symmetry properties of z/T 2 , only odd 
moments of the structure function enters the expansion of the invariant amplitude. 

Comparing equations |1.89j and |1.71| we can relate the operator product expansion 
of the forward amplitude with the moments of the (physical) structure function and 
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we get 

Au 2 E p V (p\\6f\\p) Ef(Q 2 ) = - f X dx x n uW 2 (x,Q 2 ) , n even. (1.90) 

7T JO 



That is, although the structure function 1/W2 itself can not be related directly to an 
expansion of local operators and their corresponding Wilson coefficients, its moments 
in %Bj can indeed be expressed as a product of local operators and Wilson coefficients 
by application of operator product expansion and dispersion relations. The result, 
after a trivial rewriting, is 

[ 1 dxx n uW 2 (x,Q 2 ) = (pWO^Wp) E%\Q 2 ) , neven . (1.91) 

Jo 2 . 

The Operators 

Before we discuss the momentum evolution and renormalization group properties of 
the operators and their corresponding Wilson coefficients in equation |1.91| , let us 
identify the operators themselves. 



By recalling equation |1.60| we know that the dominant operators should have the 



smallest negative dimension since those with larger negative ones will be accompanied 
by extra suppressing factors of ^ at high energy. In QCD, the leading twist operators 
are the quark operators 0* „ and the gluon operators „ as the following 

= 2 Qn^D^-D^qf 

f)G 9 pa T) T) fP 

In the quark operator, is the covariant derivative and D p = — igA^ with g 
the QCD coupling = YliA i „\ l /2 where A* is the fundamental representation of 
the SU(3) color group and A 1 ^ the gauge field. / is the flavor label and we can form 
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singlets and octets out of the flavor indices. The Dirac indices are all symmetrized. 
The dimension of the quark operator is [0*\ = —n — 2. 

The dimension of the gluon operator is also [O^] = —n — 2. The Dirac indices are 
also symmetrized. In detail, is 

Ol..., n = ~2 £ ^ 1 ^(^) <lfe (^) <A ."(^) 4 ^ < ^^&, (1-93) 

il,...,in-l 

where the z's are color indices and = <9 M + J2i 9fiijA l „ with fuj the adjoint repre- 
sentation of the 577(3) color group. 



We will not give an explicit proof of equation |1.92j , but rather present some 
motivations and explanations. 

0? ^ n Because the Dirac indices are symmetrized, only one gamma matrix can 
be used since two of them will give terms. The other indices have to be 
covariant derivatives not only to preserve gauge symmetry but also to have the 
most efficient way to get indices except for gamma matrices. For example, in 
the case of n = 2, we can have g/7^ x QfQfJ^ Qfi which has dimension of —6, or 
qfj^D^qf, which has a dimension of —4. The Wilson coefficient (E function) 
of the former will have an extra power of and thus the latter dominates. 

^Vi un The gluon field has two indices and dimension of [F] = —2, however, 
they are anti- symmetric. To symmetrize the //j indices we have to contract one 
of the two indices in each F^ u . Thus they are not as efficient as the (gluonic) 
covariant derivative T>^s. 
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1.4 Renormalization Group Analysis 

Recall that the operators need renormalization and they obey renormalization group 
equations |1.46| and |1.47| . Their corresponding Wilson coefficients obey a related equa- 



tion |1.51| so that their products, which are the n moments of the physical structure 



function, do not depend on the renormalization scale. 

1.4.1 Equations and Solutions 

We define the Parton distribution functions xPf(x,Q 2 )s by 

f'dxx-xpf^Q 2 ) = ( l dxx n m 2 {x,Q 2 ) = (pwo^wp) . 

j J J Zi ^ 

(1.94) 

At this point xP-^s seem more like mathematical objects rather than physical ones. 
However, in the Parton Model of DIS, xP^ are indeed the (momentum fraction) 
distribution functions of quarks of flavor / inside the proton Jll| . 



In leading order of the QCD coupling a(Q 2 ) (when Q 2 is large), only quark oper- 
ators come in because photons only couple to quarks directly. Thus in equation |1.94 



the O's are quark operators O* from equation 1.92 . We have chosen the normalization 
of such that 

(27rfE p (p\dl i ^Jp) = 2p, 1 ^ n (1.95) 

for quarks of flavor / in free field theory. On the other hand, in QCD O* requires 
renormalization and so a scale fi has to be introduced. There will in general be compli- 
cated Q 2 and /x 2 dependence in the reduced matrix elements and Wilson coefficients. 
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To be more specific, if fJ?/Q 2 << 1 then E£ would have terms involving /i 2 as 

El{Q\ /i 2 ) = £((oV) \ogQ 2 /^) n ) = C f + 0((a(/i 2 ) logQ 2 //?)") . (1.96) 

If we simply choose fi 2 = Q 2 , then all the Q 2 dependence is in the reduced matrix 



elements obtained by equation |1.60| and we have 



E{{Q 2 ) = e 2 + 0(a(Q 2 )) (1.97) 

under the normalization convention we use. Therefore we have, in leading order, at 
the renormalization scale fi 2 = Q 2 , 

X> 2 J^dxx^^Q 2 ) = { ^^j:( P \\df n \\p)e}, (1.98) 

which means that up to kinematic factors the moments of the parton distribution 
functions are the reduced matrix elements in operator product expansion, or explicitly, 

f Q dxx n xPf{x,Q 2 ) = ^l3L{p\\6l\\p) Q2 . (1.99) 

The subscript Q 2 means that the reduced matrix elements are renormalized at the 
scale [i 2 = Q 2 . 

Recall that the operators obey the renormalization group equation (see equations 
1.46 and 1.47] ), at an arbitrary scale /x, 



+ ^ ] d " (/i2 ' ^ = 7 ' K) ^ (/i2 ' ^ ' (L1 ° 0) 
where 7^ is the anomalous dimension of the operator. 

The matrix elements of the operators will obey the same equation. Furthermore, be- 
cause the differences between the matrix elements and the reduced matrix elements 
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are only kinematic factors (see equation |1.60| ), the reduced matrix elements and thus 
the n-moments of the parton distribution functions also obey the same renormaliza- 
tion equation, namely, we have 

(/i 2 — + (3t— ) f dxx n xP i {x : Q 2 )= 1 {{a tl ) dxx n xP f (x,Q 2 ) . (1.101) 
a/i z acta Jo Jo 



The solution to the renormalization group equation |1.100| is (formally when we have 



operator mixing and an anomalous dimension matrix) 

Qj(Ai a ,a M ) = c- J -M da ^O{(^,a, ), (1.102) 

and the moments of parton distribution functions will obey exactly the same evolution 
in momentum scale. 

The anomalous dimensions are calculable in perturbation theory. A lowest order 
calculation gives the dominant large /i 2 dependence. 



1.4.2 Relationship to DGLAP Evolution 

The DGLAP equations (Dokthitze, Gribov, Lipatov, Altarelli, Parisi) (for example, 
see state that the momentum scale Q 2 evolution of the parton distribution 

functions Pf(x, Q 2 ) obeys 

Q 2 4yi Q 2 ) = ^flf d 4 x'pf(x', Q 2 ) , (1.103) 

dQ z 2-n J x x' x' 

where r y-'(x) is the so-called Altarelli- Parisi splitting function of the corresponding 
parton distribution(s). ^{x) is defined to be zero outside the range (0, 1). Again in 
general y (x) can be a matrix due to mixing among the parton distribution functions. 
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The DGLAP evolution equations and the renormalization group equations |1. 101 



are actually equivalent, with the proper identification of the splitting function with 
the anomalous dimension. We will not explicitly prove this (for that see, for example, 
0, |[]), but rather show that given the DGLAP equation |1.103| we can obtain equation 
1 . 10 1| and illustrate the explicit relationship between the splitting function and the 



anomalous dimension (in the leading order). 

Taking the n + 1th moment (here we adopt the convention that nth moment is 



take by convoluting with x n 1 ) of both sides of equation |1.103| we obtain 

Q 2 4^2 r dx Q 2 ) = ^ t dx x n I' ^V(^yp/(a/, Q 2 ) . (1.104) 

dQ z JO Z7T JO Jx x x 

Change the integration variables in the right-hand-side (rhs) from x and x' to x' and 

x/x' by dx = x' d(x/x') and dxdx' jx' = d(x/x') dx' , so that 

£ d *l 1 .iF=C* l SXi>)< (1 - 105) 

where we have also changed the integration limits appropriately (see figure ^ for 
illustration). Rewrite x n as (x/x') n (x') n we have 

rhs = ^ f dx> f d{-ypf{x\Q 2 ){x>n-Tl f (^) 

Z7T JO JO X X X 

= ^fl [ l dxWpf{x>,Q>) / 1 rf(^)(^)V(^) 
2n jo jo x' x' x' 

_ a{gj ri dxx n xP f^ Q 2^f (11Q6) 

Z7T JO 

where we have defined the n-moment of the splitting function as 



7 ,{ = / dxx n -y f (x) . (1.107) 



By the proper identification between the n(+l)th moment of the splitting function 
and the anomalous dimension in leading order as 
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and noting that at the normalization scale /i 2 = Q 2 we have a M = a(Q 2 ) and thus 

= (1 ' 109) 



we have indeed recovered equation |1 . 101 



1.4.3 The Anomalous Dimensions and Splitting Functions 

In terms of diagrams, the operator product expansion of the DIS forward amplitude 
into the product of reduced matrix elements and Wilson coefficients is shown in 
figure [7]. If we put a cut through both sides, it becomes the OPE of the structure 
function into Wilson coefficient and parton distribution function. Only the lowest 
order contribution is shown where the connecting part on the right hand side is a 
quark operator. 

The parton distribution on the bottom evolves in the momentum scale via an anoma- 
lous dimension until it reaches the scale of the hard scattering Q 2 . We will continue 
the computation and discussion in the light cone (LC) gauge. 

Review of the Light Cone Gauge 

We take the convention that the component is the time component of a four vector, 
and the 1,2,3 components are the (spacial) x,y,z components, respectively. 

For any four vector v, we define the light cone components 



1 . 

v± = y|( u o ±1%), v = 



U 

\ v 2 



1.110) 
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This means 

v 2 = v ■ v = 2v + V- — v 2 , (1.111) 
while for any two 4- vectors V\ and v 2 , 

v 1 -v 2 — v 1+ v 2 - + v l _v 2+ - 2h ■ v 2 . (1.112) 

We can then define the LC metric g a p by 

9+- = 9-+ = 1, 9++ = 9— = 0, 9±± = -1, 

0+± = #-L+ = 0-± = = 0. 



such that for a, (5 — +, — , _L 



«1 • «2 = J^ffa/S^la^/?- (1.113) 

a,/3 



7± = ^(7o±73), 7 = 



1.114) 



The Dirac matrices 7 a with ct = 0,1,2,3 can be viewed as a four vector and we 
can similarly define 

/ 

7i 

V 72 

It can be easily verified that the light cone components obey the anti-commutation 
relationship 

{7a, 7/3} = 2g af} . (1.115) 



In particular, we have 



7+ = 7! = 

7+7- + 7-7+ = 2 (1.116) 
7±7i,2 + 7i,27± = 



and therefore, 

7+7-7+ = 27+ 
7+ il+ = 2^+7+ 

Introducing a vector such that for any 4-vector v 

n • v — v + . 

This means 

n_ = 1 , n + = n± = , & n 2 = n^rf" = . 

n M is the so called LC null vector. 

The light cone (LC) gauge is defined by requiring for the gauge field 

n ■ A = A + = . 
The gluon propagator in the light cone gauge is 

Ar + «e \ n ■ k J k z + ie 

where D^ v (k) is the so-called light cone gluon projector. We have 

n^ u = n u V^ = . 

The light cone gauge is also a ghostless gauge. 
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The Quark-Quark Anomalous Dimension 

The diagram of the quark distribution function xP*{x) is shown in figure || where 
we take the lowest order quark-quark vertex function to be ( ||45|| ) 

rf = 7+ g±)". (i.i23) 

Let us look at the first radiative corrections. The quark-quark diagram is shown 
in figure |9], where we have an extra on-shell gluon line. 

Using standard Feynman rules we can write the value of part of the diagram within 
the dotted circle as 

In the above, C is the color factor and 

C = (^) 2 E(^aO, (1-125) 

i.b 

where T % ah = ^ is the generator of the SU(3) color group. Since we have 

]T(r) 2 = C F /, (1.126) 

i 

where CV is the Casmir operator of the fundamental representation of the gauge 
group and for SU(N) 

N 2 - 1 

C F = - -. (1.127) 

2N K J 

Defining the strong coupling constant a s = a by g 2 = Aira, we have 

C = (ig) 2 C F 5 a , a = -AnaCp , (1. 128) 
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where in the last step we have set a = a'. 



In arriving at equation 1.124 we have also used the transition from the usual gluon 



propagator |1.121| to the on-shell gluon by 



k 2 + ie 



(-27r)5+[fc 21 



1.129) 



which is consistent with the optical theorem |1 .39 . 

Therefore we obtain 

(1) _ aC F f d A k5 + [(k 1 - A;) 2 ] (k 



n— 1 



2tt 2 



7o 



^D a/3 (h-k). (1.130) 



(k 2 + ie) 2 \p 4 

Note in the light cone gauge, there is only one solution to 5(k 2 ) because of the 
linearization of the light cone variable. We write the integral over d A k in light cone 
variables 

d 4 k = dkodkidk 2 dk 3 = dk + dk-d 2 k, (1.131) 



and integrate over k- with the 5-function to get 



dk^Kh-k) 2 ] = dk.sfiih-ty+ih-ty^-^-kf] 



2(k 1 -k). 



1.132) 



Define three momentum ratio variables 



k, 



x 



Xi 



k u 



x 



we evaluate k 2 with {k\ — k) 2 = and get 



ki + x\ 



1.133) 



k 



2k + k_ - k 2 
2k + (k v 



{k.-kf 



2(k 1 -k). 



kt^ + k ^ k * (h-k)^ 1 - )2 ' 



1.134) 
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where in the last step we have used 

k-u . k 



2k+h_ = ^2k 1+ h_ = P-{k\ + fc?) . (1.135) 
h+ k l+ 



Expressing things in terms of u we have 



k 2 = uk 2 1 +uk 2 -k 2 --^(k 1 -k) 2 

1—UJ 



1 [-u(l-u)k% + uj^-k) 2 - uj(l-u)kl + (l-u)k 2 



= --^—[{k-ujk-tf-uj^-^k 2 }. (1.136) 

1 — LU 

Our goal is not to calculate T + completely, but only to calculate the ultraviolet 
divergent part of it. In evaluating the numerator of |1.130| we keep only the terms of 
the highest k power, which turns out to be quadratic. Note because of the 5-function, 
(&4 — k)~ also has a A; 2 contribution. The denominator is simple after we take the 
leading contribution, that is, 

k 2 = --^-. (1.137) 

1—UJ 

As for the numerator, we now need 



7^7+^7^ = 7 Q ({^,7 + }-7 + W VrfDafi 

= 2k +1 a ft-fDvp - k 2 ^ a j + ^D a/3 . (1.138) 



The second term of the above is 



n ■ (ki — k) 



k 2 

2k 2 1+ = -2 7+ -=- , (1.139) 

1 — LU 



47 



where in the first step we used the facts that n ■ 7 = 7+ and 7+ = while in the last 



step we used the anti-commuting relation |1.115| . The first term of the numerator, on 
the other hand, can be expanded to be 

2k + { la #7° - * [7+ Mi- W + ih~ W ■ (1-140) 



We have from the 5-function 



k -~ ' (L141) 



and thus 

k 2 



la ^ a = -2f6= -2 1+ k_ = 7+ ~ . (1.142) 

{ki-k) + 



At the same time, 

7+ = -l + ^=~l + k 

k 2 

l-UJ 

while similarly 



7+A~ . (1-143) 



fx- ft #7+ = 7 +T ^ ■ (1-144) 

1 — C<J 



Therefore the first term of the numerator is actually equal to 



2 U , 2 2 UJ 1+UJ 



27k 2 - {1 - } = -2-yk*- . (1.145) 



Combined with the second term, we find that the divergent contribution to the nu- 
merator is 

k 2 r . , . „ „ , 2 l+U 2 



~ 2l+ T^uf [uj{l+uj) + = " 27+ - (WP ' (1 ' 146) 

Rewrite the integral as 

d 2 £ = 2vr / kdk = n f dk 2 , (1.147) 
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we arrive at the expression of the divergent part of r+ as 

= i „ n dk \, {-V i^iA^-Vi— y 1 2(1+ up 



2tt J 2(k x - k)+ \p + J J k 2 \k 1+ J \p+ 

,'ki+\ n ~ aCp [ dk 2 r , „ A+uj 2 

7+ — / -TT / duuo n ~ 1 ^— . 1.148 

1 1 p + J 2tt J k 2 J 1-uj 



Our original graph, one that is figure 0(b) plus an on-shell gluon running from 
the incoming quark to the outgoing quark, is not divergent in k 2 . The divergence we 
encountered arises because we let Q 2 become very large for fixed fc 2 . This means that 
the cutoff for the logarithmic fc 2 integral is at Q 2 . That is, explicitly, 



where we have also explicitly written out the limit on the u integration. 

Taking the logarithmic derivative with respect to Q 2 , and recall equation |1 . 123 
we obviously have 

Q 2 -L^M 2 ) = ^ Cduu^-^rf (x,Q 2 ). (1.150) 



dQ 2 + v ' 2tt Jo l-u 

Remember the graphic definition of the parton distribution function and as 
the lowest and first order factors, we realize that we indeed have 

Q 2 4t2 t^xPfix^Q 2 ) = ^ /V^P'foQ 2 ) tduu^ 1 -^. (1.151) 

dQ z JO Z7T JO JO 1—LO 



By exactly reversing the argument in section |1.4.2| (with x' now being xi), it is straight 



forward to show that we indeed have the DGLAP evolution equation for the quark 
distribution function 

Q 2 ^xP f (x, Q 2 ) = ^ L [ 1 ^7(-)xP / (x 1 , Q 2 ) , (1.152) 
aQ z Zir J x x\ x\ 
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where the Altarelli-Parisi splitting function is 



l + v 2 

i(y) = ^ L . (1.153) 
1 -y 



The integral in equation |1.152| is actually divergent at the end point when x\ — > x. 
This is an infrared divergence due to emission of soft gluons with {k\ — k) + — > 0. This 
divergence should not be present in physical processes and indeed it is removed once 
we add the quark self-energy graphs shown in figure |IJj. 

After virtual corrections are included, the full quark-quark splitting function that is 
free of singularities is 

j(x) = C F [^±^- + ^8(l-x)) 1 (1.154) 
where the special function (1 — x)^ 1 is defined by 

/(,)* _ fyw-/(i) (L155) 



/o (1 — x) + Jo (1 — x) 
for any function f(x) that has reasonable behavior. 

The Mixing of Evolution 

There are additional order a terms in the Q 2 -evolution equation due to operator 



mixing, for example, see figure [IT]. We thus have the necessity of considering both 
quark and gluon distributions. 

To organize, define the flavor singlet quark distribution as 

S(a;,Q 2 ) = xJ2[P f (x,Q 2 ) + P f (x,Q 2 )], (1.156) 
/ 
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while the flavor octet distribution is 



A f f'(x, Q 2 ) = x[P f {x, Q 2 ) - P f '{x, Q 2 )} 



(1.157) 



and a similarly defined function A^' for anti-quarks. 

The octet distributions obey a generic DGLAP evolution equation 



(1.158) 



However, the singlet distribution will mix under renormalization with the gluon dis- 
tribution G(x,Q 2 ). The DGLAP equation has a matrix form and explicitly 



d I S(x,Q 2 ) \ _a(Q 2 ) r l dxi I Iqq'^/xx) -? qg (x/xi) 
d Q 2 \G(x,Q 2 )) 27T ■** x ^ \i gq {x/ Xl ) lgg {x/ Xl ) 



S(x 1; Q 2 ) 
G(x^Q 2 ) 



1.159) 



Let Cf again be the Casmir operator in the fundamental representation with Cp = 
4/3 for SU(3) color group, and let Ca be the Casmir operator in the adjoint repre- 
sentation with Ca = N c = 3 for the color group, we have the explicit expression of 
the four splitting functions as (eg, see [0]) 



lgg( x ) 



C F 

C F 

1 r 
2 



l + x 2 3 , 
^xY + + 2^- x) 
l + (l-x) 2 



x 2 + (l-x) 2 



2C A 



X 



1 — x \ 
+ + x(l-x) 

(1— X) + X 



(1.160) 

(1.161) 
(1.162) 

UC A -2n f 

J -8(l-x) (1.163) 

6 
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1.5 High Energy Behavior 

It is more convenient to discuss the high energy behavior of the parton distribution 
functions when we go to the n-moment space by taking the n-th moment of the 
DGLAP evolution equation, or equivalently directly using the renormalization group 
equation from operator product expansion. 



1.5.1 Anomalous Dimension Matrix 



To the leading order in a we obtain, after taking the moments of the splitting func- 
tions, the anomalous dimension matrix for the singlet-gluon mixing to be 



ln(a(Q 2 )) = * ) 

fn 9 ) 



2tt 



1.164) 



in \r< 



where 7„ is the n-th moment of the corresponding ^y(x) in equation |1.16(J| and explicitly 

2 



7 



ry19 
I n 



ry99 
I ii 



ry99 



n i 

1 + 4V-- 

2 3 n \ n + 1) 

n 2 + n + 2 



n(n + l)(n + 2) 
In 2 + n + 2 



n 


n 2 - 


-1) 




"1 


2 


3 


3 + 


9^ 



1 



^ j n(n-l) (n + l)(n + 2) 
1.5.2 Dominant Moment Contribution 



1.165) 
1.166) 
1.167) 
1.168) 



The anomalous dimension matrix elements are usually singular at some values of 
the moment index n (see equation 1.165 ). In the high energy limit, the dominant 
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contribution at small x is from the right most pole of the anomalous dimension. We 
will illustrate this by looking again at the DIS structure function uW 2 . 

Adopting the common terminology let us now call vW 2 structure function F 2 . The 
n-th moment of F 2 is defined as 

F 2 (n) (Q 2 ) = f 1 dxx n ~ 1 F 2 (x, Q 2 ) . (1.169) 
Jo 

As we already established can be factorized as the product of parton distribu- 
tion functions and perturbatively calculable Wilson coefficients with well defined Q 2 
evolution. To predict cross-section, however, we do need directly F 2 (x, Q 2 ) itself. We 
can reconstruct F 2 from its moments by a Mellin transformation 

F2(x,Q 2 )= f^.e n ^Ft\Q 2 ), (1.170) 

Jc 2,711 

where c is a contour lies in the complex n plane parallel to the imaginary axis and to 
the right of all singularities. 

Let us first check to see whether this is a self-consistent definition. Substitute 



1.169| into equation |1 . 1 70| we obtain for the right hand side (rhs) 

dn r 1 



rhs = { [ dx'x'^F^x'.Q^e^ 1 - 
Jc 2rn Jo 

= I dn f 1 dx^ e nlo g l-nlo g j TF ^ x ^ Q 2 ) 

Jc 2"7U Jo x' 

= f 1 ^F 2 (x', Q 2 ) ( ^e^o^-io^] 
Jo x' Jc 2wi 

= f 1 d\ogx'F 2 (x',Q 2 )6(\og--\og^ 

Jo V x x j 

= F 2 (x,Q 2 ) = lhs, (1.171) 
where in the fourth step we have used the mathematical identity 

L—ioo 27T2 
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We can always distort the contour c to the left (but not crossing any poles) so that 
it'll pick out the residues of the n poles. Each pole n p will eventually translate to a 
power factor of the Bjorken-x as (x)~ np . Therefore, at high energy and small- £ the 
dominant contribution comes from the right most n pole of the anomalous dimension 
on the complex n-plane. 



1.5.3 Double Leading Logarithmic Approximation 

To illustrate the point that leading (right most) pole in moment space of the anoma- 
lous dimension dominates in the high energy limit, we will explicitly compute the 
high energy behavior of the structure function F 2 of deeply inelastic scattering. 

According to previous analysis, the moments of the structure function obey the 
momentum evolution equation (see, eg, |1.102| ) 



Ft\Q 2 ) = Ft\Ql)J*l ^ 7 " WA2)) . (1.173) 



Thus the structure function itself is given by 



2 



F,(x,<?) = I >(Q3) e **i + £8 . (1.174) 



From equation |1.165| it is obvious that the gluon-gluon anomalous dimension has the 



right most pole among others at n — 1. Using the lowest order running coupling we 
can write the leading contributing term of the anomalous dimension 

TnHA 2 )) = , Ca Uh = \ 2/A2 , (1.175) 
n(n — 1) b log A 2 / A 2 

where 6q is the first coefficient of the QCD /3-function. 
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By substitution we have 



F 2 (x,Q 2 



I 
I 



dn 
2ni 



Ft\Q 2 )e" 



x 7t6q ° i„„ cfl /a2 



logQ2 /A 2 n -l 



(1.176) 



where we have defined at high energy, 



a = log — > > 1 

x 

Ca logQVA 2 

" TTbo 0g l0gQg/A 2 



(1.177) 



The high energy limit corresponds to n — > 1 and a » 1 so we can make a saddle 
point approximation for the exponent function f{n) = an + &^ty- We have 



f(n) = a- 



Thus at the saddle point n , 



(n-l)* 



=^ n = 1 



(1.178) 



- + 1) + -= = + 2^ 
/"(no) = 



2a 3 / 2 

6 l/2 



>0, 



(1.179) 



which means c is a fine contour. Let iv = n — 1 we obtain finally 



F 2 (*,Q 2 ) = Ft\Ql) T d ^^e-t 
= F 2 {n \Ql)e n ^e 



3/2 

v 2 - 



V a 3 / 2 



= F 2 W (Q 2 ) ( J ) e 



\ (log7) 3/2 



(1.180) 
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This is the high energy limit of the forward DIS amplitude/structure function. 
The dominate small- a; behavior is a power dependence on ( with n given by a 
saddle point approximation. It is obvious that the value of n , which determines the 
leading high energy and small-a; behavior, is dictated by the right most pole on the 
complex n plane. In this forward case, it is at n — 1 and thus we arrive at the x~ l 
leading behavior of the gluon distribution under DLLA. 
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Chapter 2 

Decomposition of the Non-forward 
Amplitude 



In this chapter, we discuss in detail the tensorial decomposition of the non-forward 
Amplitude. First we define the amplitude under discussion by laying out the kine- 
matics of the process; then we proceed to use current conservation and the symmetry 
properties of the amplitude itself to decompose it into various invariant components 
and define the corresponding invariant amplitudes in the non-forward case. 

2.1 Kinematics 

The process we consider is 7* + P — > 7* + P (virtual photon + proton goes to virtual 



photon + proton) shown in fig. [12] with its amplitude. It is a double virtual 
Compton scattering in the sense that both the incoming photon q' and the outgoing 
photon q have, in general, non-vanishing invariant masses. There is a non-zero four- 
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momentum transfer from the virtual photon to the target proton which we label as r. 
Thus the struck proton p' has a slightly different four-momentum than the original 
proton p. And we have 

p = p — r ; q = q — r . (2.1) 

Furthermore, the incoming virtual photon q' has a different invariant mass than the 
outgoing virtual photon q with 

q 2 = ~Ql q 2 = ~Ql (2.2) 

We restrict our discussion to amplitudes. Cross sections are obtained from the 
square of the amplitudes. In this non-forward case there is not a simple optical 
theorem relating the imaginary part of an amplitude to a cross section |12j| . Our 
main interest is in the high energy limit, so similar to [|l7| we take all but one of the 
light cone (LC) components of the proton momentum to be zero. We do the same 
for the momentum transfer r. We will choose the non-zero component of both as the 
plus component, which means 

p=(p+, 0,0), r=(r+, 0,0), (2.3) 

and 



p 2 = = r 2 = -t . (2.4) 

That is, we are in the zero nucleon mass and zero t limit. 
The fact that the target remains a proton means that 

p' 2 = = {p — r) 2 = p 2 — 2p ■ r + r 2 , (2.5) 
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and thus 



which gives 



And since 



which results in 



p-r = 0, (2.6) 



p. q ' = p.q. (2.7) 



/2 / \ 2 2 22 /\ 

g = (g — r) =g — 2g • r + r =g — 2q ■ r , (2.8) 



2g ■ r = q 2 - q' 2 . (2.9) 



Together with the fact that r is proportional to the external momentum p, we 
have, 

r = (p, C = 4^- (2-10) 
2p ■ q 

We always suppose — q 2 — Q\<Q\ — q' 2 so that ( > 0. 

It is worth noting that the process of a double virtual Compton scattering is 
not a physical one, however, it provides a general framework in which one can move 
the discussion continuously from one physical limit to another. In particular, DIS 
corresponds to q — q' while DVCS corresponds to q — 0. 

Also worth noting is that the kinematic limit we are discussing, namely, r ^ but 
t — 0, is not a physical one, since for a physical non-forward process, t is bounded by 

x 2 m 2 

\t\>—, (2.11) 

where x is the Bjorken x variable and m the target proton mass. However, our main 
interest lies in the high energy and small-x behavior of the amplitude and therefore 
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t = is a good approximation. In principle we can always discuss finite t behavior 
by giving r a transverse component r ^ 0. There the proportionality between r 
and p (see [2.10[ ) becomes a proportionality between the plus components of the two 
momenta (i.e. r + = (p + )(see, for example, f23fl). 



2.2 The Amplitude 



The amplitude of the process in fig. [12] in coordinate space is a matrix element between 
two proton states of a time ordered product of two electromagnetic currents. If we 
label the space time point of the incoming interaction as x + z and the outgoing one 
as x, we have 

T coor d = tf\ Tjii ( x ) jv ( x + z )\ p ) . (2.12) 

As usual we want to Fourier transform the amplitude into the momentum space where 
the kinematics and the dynamics of the process can be much more readily discussed. 
The momentum space amplitude, which we label as T^, depends on all the relevant 
independent kinematic variables (in the momentum space). In this non-forward case 
there are three kinematic degrees of freedom and sets of variables such as (p,p', q) or 
(p, q, r) or (p, q', r) are valid choices and are all equivalent to each other. We choose 
the set (p, q', q) for convenience of later discussion. Therefore, using the convention 
that an incoming momentum p at space time point x enters the Fourier transformation 
as an exponential of negative exponent e~ tp ' x and an outgoing momentum as one with 



a positive exponent, we can write explicitly the amplitude of fig. |T| in momentum 

space as 

TV„ = T^(p,q',q) 
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% f d 4 xd 4 z e -V ( p '\T j„{x)j v {x + z)\p) . (2.13) 



Integration over d 4 x gives only an overall four dimensional momentum space 5- 
function. Ignoring it we can rewrite the amplitude, up to normalization factors like 
(2tt) 4 , as 

T^(p,q',q) = i J d A ze^'- z (p'\T j^)j v {z)\p) . (2.14) 

It is obvious that this amplitude should be invariant under an overall translation in 
coordinate space, for example, (z — > & — > —z). Explicitly, by using the translation 
operator e l ^' x where p is the momentum operator and for any operator in coordinate 



space 0(z) (see equations |1.18| and |1.19|) 



e l ?- x d(z - x)e~ ip - x = 6{z) , (2.15) 



we have 



TrM,q) = tJd*ze-^ z (p'\Tj,(0)Uz)\p) 

= i [ d 4 ze~ iq '- z (p'\Te ip - z j^-z)e~ ip - z e l P- z j u (0)e-^ z \p) 



-ip-z 



i J d 4 ze-^- z e^- z (p'\Tj,(-z)U0)\p)e- 
i J d 4 ze-««'-^>*(p'\T U-z)j v (0)\p) 
i [ d 4 ze-^ z (p'\TU-z)U0)\p) . (2.16) 



Anticipating the operator product expansion analysis we want to recast the amplitude 
into the so-called light cone expansion where [z — > ~z & — > —\z). Similar to the 
above process we have 

TV(p,</,g) = % fd 4 ze^'- z (p'\TU0)Uz)\p} 
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= i J d 4 ze- ig '- z (p'\T e ip ^ z j^-^)e- l P^ z e ip ^ z j u (^z)e- if, ^ z \p) 

= ijd'ze-M+^ip'lTj^-z^mp). (2.17) 

Therefore, we arrive at 

T,u(p,q',q) = i J d'ze-^ z tf\T ^(-|)>(|)b) (2-18) 

where we have defined 

g = g' + -r = g--r = -(g' + g). (2.19) 

Because from ( |2.9| ) we have 

g ■ g' = g 2 - g • r = i(g' 2 + g 2 ) = -1(Q 2 + Q 2 ) , (2.20) 

the square of this newly defined g becomes 

q 2 = \(q' + q) 2 = \(q' 2 + q 2 ) = q-q'. (2-21) 

If we define 

Q 2 = ~t = -\{q' 2 + q 2 ) 

= \{Ql + Ql), (2-22) 



we will see in later discussions that Q is now the natural scale of the scattering 
process and it characterizes the hardness of the scattering. 
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2.3 Current Conservation 

Similar to the case shown in the previous chapter for the forward scattering we will use 
the conservation of the electrcnnagnetic(e&m) current ( |2.23| ) to simplify the tensorial 



decomposition of the non- forward amplitude T^. However, because now the two e&m 
currents carries different momenta, we need to be more careful in writing down the 
consequences of the current conservation in momentum space. 



Because 



d Uz) = (2.23) 



from ( |2.16|) we have 

= ifdtze-V-'WTjrW—jMlp) 
= iJ*ze^'*^V\TjMM\p) 



-i J d A z{-iq') v e^- z (jp'\T j,(0)Uz)\p)) 

+ iq'jT^v (2.24) 



and 



i fd i ze-^(p'\T(^U-z))UO)\p) 
l Jd 4 ze-^—(p'\T 3ll (-z)UO)\p) 

ijfzj^le-*' ( P '\Tj„(0)j u (z)\p)) 

-i J d A z(-iq) v e-^( P '\T 3 ,(-z) ]v mv)) 

+ iq^T^ . (2.25) 
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Thus in momentum space the conservation of the electrcunagnetic current now re- 
quires the amplitude to satisfy (cf |1.22| ) 

g^V = 0; = 0. (2.26) 



2.4 Decomposition and Invariant Amplitudes 

To obtain the complete tensorial decomposition of T^ we also need the symmetry 
properties of the amplitude. From fig. |X3] we can see clearly that the second diagram 
is simply the first diagram with a different momentum labeling and thus the two have 
identical values, both equal to T^ u . Note that to be completely explicit, each diagram 



in fig. [13] should have its own corresponding cross diagram, and the symmetry exists 
for the sum of the cross diagram and the original one. 

Therefore, while T^ is no longer symmetric, it is invariant under the transfor- 
mation <-> z/, q <-> — q', that is, 

V(p,g',g) = T^(p,-g,-g / ). (2.27) 



From Section |2.2| we know that T Mi , depends on the momenta p, q and q' } similar 
to the last chapter we can write the most general tensorial decomposition of T MJ/ as 

TV„(p, q', q) = a g^ + aip^p v + a 2 q' pL q' u + a 3 q^q u + a 4 p^q' u + a 5 q^p u 

+a 6 p IM qv + a-jq^pv + a 8 q'^q u + a 9 q^q' u . (2.28) 



As in the previous chapter, this is a spin averaged electron-scattering, thus we do 
not need to consider any helicity issues and hence there are no e-tensor terms in the 
decomposition. 
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The (ZjS are invariant amplitudes depending again only on scalar products of the 
momenta. Because we have equations Q2.6|), ( |2.7| ), ( |2.9| ), and ( |2.20| ), we take the ajS 
to only depend on q' 2 , q 2 and p ■ q, that is 

a i = a i (q' 2 ,q 2 , P -q), i = 0,1, 2,..., 9. (2.29) 



From (|2T27|) and ( ^28|) we have 



TV M (P, -?') = ao{~)g^ + ai(-)p^ Pv + a 2 (-)q'^q u + a 3 (-)g M ^ - a 4 (-)p^ql 
-a 5 (-)g^ - a&(-)Pn<lv - ai{-)%Pv + as(-)^9i/ + ag(-)?X 
= TVfog'.ff) (2.30) 

where we have used the short handed notation 

Oi(-) = -</) = a;(</ 2 , g 2 , -p-q) ■ (2.31) 



Equation |2.30| is true for any values of the momenta g', g and p, in particular, we can 



fix the scalar products but still vary arbitrarily the individual momentum. Thus the 
coefficients of the same momentum combinations must be identical, which leads to 
the following relationships among the invariant amplitudes: 

ao = ao(-) ; «i = oi(-) ; a 8 = a 8 (-) ; a 9 = a 9 (-) ; 
«2 = «3(-) ; a 3 = a 2 (-) ; a 4 = -cn(-) ] «5 = -ae(-) • ( 2 -32) 

By using the conservation of the e&m current ( j2.26j ) we can establish another 
two sets of equations of the invariants. We have 

q^T^ = = a q u + a x p ■ qp v + a 2 q ■ qq v + a 3 q 2 q u + a 4 p ■ qq v + a 5 p ■ qp v 



65 



+a 6 p ■ qq u + a 7 q 2 p u + a 8 q ■ qq v + a 9 q 2 q' u ; 



q' u T^u = = aoq'^ + a\p ■ q'p^ + a 2 q q^ + a 3 q ■ q'q^ + a 4 q p^ + a 5 p ■ qq^ 



+a 6 <? ■ q'Pn + a 7 p ■ 4% + a s q ■ q'q' + a 9 q' g M . 



(2.33) 



Because again this is always true when we vary arbitrarily the momenta while keeping 
the scalar products fixed, the coefficients of the same combinations of momenta must 
identically vanish, leading to 



Oo + a 3 g + a 6 p- q + a s q ■ q 
a 7 q 2 + aip- q + a 5 q ■ q 
a 9 q 2 + a 4 p ■ q + a 2 q ■ q 
a + a 2 q' 2 + a 5 p ■ q' + a 8 q ■ q 
a 4 q' 2 + dip ■ q' + a 6 q ■ q' 
a g q' 2 + a 7 p ■ q' + a 3 q- q' 












(2.34) 



As in the previous chapter, we will use these equations to reduce the number of 
independent invariant amplitudes. We first rewrite ( |2.34j ) by taking 00,1,2,3,5 as given 
and expressing the rest of a'^s in terms of them. We obtain 

,2 



a>4 



q • q 

p-q 

p-q ' 

r ,l2 



-a 2 



q 



-ai — 



p-q 

q ■ q' 



-a g 



a 6 



p-q 



a 7 



q ■ q' 

p-q 
p-q 



7 ai 



p-q 

q' 2 



a 2 + 



1 1 

q ■ qq 

q' 2 p ■ q 



a 3 - 



q ■ q 



V2 



-a 5 



-a 2 



q ■ q 

q 2 



-a 4 



q z 



-Oi 



p-q 
q-q' 



a 3 + a 5 



-a 5 . 
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,'2 



I q'" p ■ q' 



a g 



q ■ q 
Q ■ <t 

■ q ,2 



-a 3 



q ■ q' 

p-q ' 

q' 2 



q ■ q' 



-a-j 



p ■ qp ■ q' q ■ q' p ■ q'q ■ q' 

-o-i ^-«3 H — a 5 



q 2 q' 2 



q 2 q' 2 



(2.35) 



We now use equation 2.32 to further simplify the above result. We have 



a 8 (- 



-/2 



g • q 



I g'" p • g' 
; a o 7 a 2 



9-9 



g • g' 



g-g 



7 a o(" 



g-g' g-g' 



p ■ q 

«5 + a 5 (-) = —^-(03 - a 2 ) 



a 9 



Og(-) 
p ■ qp ■ q' 



q 2 q' 2 



-CL\ 



q ■ q 



q 

p ■ qp ■ q 

q 2 q' 2 



r'2 



a 3 + 



p ■ qq ■ q 

2~i2 ° 5 



p ■ q'q - q' 

q 2 q' 2 



«5 + a 5 (-) = ^r(a3 - a 2 ) 
g z 



(2.36) 



where we have used a\ = a\(— ), a 3 (— ) = a 2 and p ■ q = p ■ q'. From the above, it is 
obvious that we must have 



a 2 = a 3 , a 5 = -a 5 (-). 



(2.37) 



Furthermore, 



a 5 



-ae(-) 

a 5 = -( — - — a 2 (-) + a 3 (-) + a 5 (-)) 



-,'2 



a 5 + a 5 (- 



p ■ g p ■ q 

g' 2 g 2 

a 3 a 2 

p • g p ■ q 



(2.38) 



which, combined with (|2.37|) , gives 



a 2 = a 3 = . 



(2.39) 
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At the same time, 

04 = — aj(— ) =>- 

p-q' q-q' q-q'q 2 q ■ q' p ■ q q ■ q' 

^-ai a 2 + — a 3 —a 5 = ai(-) — a 5 (-)) 

q/z p . g q/zp . q q/z qz qz 

=> a 5 = -^-l ai . (2.40) 
Therefore, going back to ( |2.35| ) we have 



p-q' q-q' p-q 
q/z q/z q . qi 







p-q 

a 6 — °5 — 7 a l j 

q-q' 

p-q q-q' p-q 

a 7 = 5~ai 5~( -Mi 

q A q q • q 







1 p-q' p-q 

a 8 = -oq -( -Jai 

q-q q-q q-q 



1 (p-q) 



2 



p-qp-q' V-q'q-q' p-q. 

a 9 = Y7z~ ai + Fa - ( 

g z g w g z g z q ■ q 

= 0. (2.41) 



It is now clear that only two independent invariant amplitudes, a and ai, remain 
in the end while all the other ones either vanish identically or can be expressed in 
terms of ao and d\. In summary, we have 

0-2 — a 3 = a 4 — 0-7 — a 9 — 

p-q 

a 5 = a e = -ai 

q-q' 

a 8 = —a + ^^ai . (2.42) 

q-q' (q- q') 2 
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Thus, substitute these back into (2.28) we have 



P ' 9 

TVO, q , q) = a g^ + a x p tl p v - ——^{q^ + p^q u ) 

,/ 1 , (p-qT s , 
9 ■ q \q ■ ? ) 

/ 9^ \ 

, ( v-q , / , x . (p-9) 2 , \ / 0/l ^ 

By defining two invariant amplitudes T\ and T 2 we finally obtain the general 
tensorial decomposition of as 

1 rp / , 9^\ 



+ ^ - ^-(P^u + q'^Pu) + (^)V„& ) T 2 . (2.44) 



T; = Tj(g /2 , q 2 ,p ■ q) are invariant amplitudes analogous to those in DIS (see section 
1.3.3|) . They are even functions of p • q (because they are essentially and a{). 



We have pulled out from Ti an explicit factor of y/1 — (, as was done in, e.g. ,|L6| (see 
equation [2.10| for the definition of () . This factor comes from the external spinors 
of the proton as the following: The normalization convention we use for computing 
Feynman diagrams is such that an external incoming Fermion with momentum p 
and spin state r enters a Feynman diagram expression as Ur ^ where U T (p) is the 
standard 4-component basis spinor for a Fermion. Thus there is an explicit factor 
of y/p7 1 = (VI — Cy/p) ~ 1 i* 1 the expression of T^. We pull it outside explicitly so 
that later on the expression of the lowest order Wilson coefficient will be simpler (see 



equation 4.60 in section fl 



It is worth noting that we still have only two independent invariant amplitudes 
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even in the presence of three independent invariants, in contrast to the case of forward 



scattering where we only have two (see sections |1.2.1| and |1.3.3| ). 
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Chapter 3 

Operator Product Expansion 



In this chapter we will discuss the operator product expansion (OPE) of the non- 



forward amplitude (see section |2.2|) . We will identify the operators in the OPE, 
obtain their reduced matrix elements between (asymmetric) external states and define 
new moment variables. We will express and eventually the invariant amplitudes 
Ti and T2 (see equation p.44|) in terms of the reduced matrix elements and the 
corresponding Wilson coefficients. 



3.1 Non- forward Operator Product Expansion 

In the short distance limit, we can perform an OPE for as a sum of products of 
local operators and their corresponding Wilson coefficients ||: 

y y n oo J Uj w 

Z Z J=ln=0i=l Z Z 

00 J+2 uj ix x uj 

+ E E E e%{z^)6%^^) z —... z — , (3.1) 

,/=l 71=0 8=1 Z Z 
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where A^ v and B^ ua p are conserved tensor structure operators. The summations 
should be regarded as such that in addition to the literal summation over J, n and i, 
all the indices except ji, v are also symmetrized as usual. 

To obtain the explicit expression of B^ a/ 3 we start with the most general form 
that an tensor operator in the position of B can have: 

+a49 f Mud a d (3 + a^g m d v dfi + a 6 g^d a d v 

+a 7 g u(3 d^d a + a 8 g ua d^d p + a 9 g a/3 d^d v . (3.2) 

Current conservation requires 

dj^O^ap = dyO^ap = . (3.3) 

We have 

= aid u g a(3 U + a 2 d a g v pU + a 3 d p g ua a + a$ v d a dp 

+a 5 d a d v d/3 + a Q dpd a d v + a 7 g u/3 d a a + a 8 g pa Ud p + a 9 g aj3 Ud v 
= ((ai + a 9 )g a pd u + (a 2 + a 7 )g u(3 d a + (a 3 + a 8 )g va dp)U 

+ (a 4 + a 5 + a 6 )d u d a d (3 , (3.4) 



and 



aid^g^D + a 2 g iia d fi U + a z g^d a U + a A d^d a dp 

+a 5 g^ a d/30 + a 6 g^d a a + a 7 d p d^d a + a 8 d a d^d p + a 9 g ap nd^ 

((ax + a 9 )g afS d^ + (a 2 + a^g^dp + (a 3 + a 6 )g^d a )a 

+(a 4 + a 7 + a 8 )d IJi d a dp. (3.5) 
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Since the above two equations are tensor equations, the coefficient of each (different) 
tensor structure operator must vanish, which gives us 

ai + a g = a 2 + a 7 = a 2 + a 5 = a 3 + a 8 = a 3 + a 6 = 

04 + 05 + as = 04 + a-j + as = . (3-6) 

Thus we have 

a 9 = -a x 
a 4 = a 2 + a 3 
0-5 = 0-7 = —0-2 

a% = a s = -a 3 , (3.7) 

which, after substitution back into (|3.2|) , leads us to 

O^ap = aig^g ap U + a 2 g m g ufi U + a^g^9va u 

+ (a 2 + a 3 )g^d a d[3 - a 2 g IMX d v dp - a z g^d a d v 
—^gvpd^da — a 3 g va d^dp — aig a pd^d v 
= {g^g a f3 u - g a f3d^d u )ai + (g m gv/3 n + g^d a dp 
-g^adudp - gupd^d a )a 2 + (g^g ua a 

+g^d a dp - g^d a d v - g^dfl^a?, . (3.8) 

If we symmetrize the indices a and (3 we have 

0^ a p = g a p(g^ u - d^d v )ai 

+{g^gvp u + g^dadp - g m d v d$ - g v pd ll d a )a 2 . (3.9) 
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It is now obvious that B^ va p is the same as in the forward case (see equation 1.54| ) 



i.e. , 

Byvap = g m g u p u + g^udadp - g m d v dp - g v pd,J) a , (3.10) 



and corresponds to the tensor structure multiplying T 2 in Q2.44Q . 



We have not found the explicit form of A^ v which would generate the tensor 
structure corresponding to TV But, as in the forward case, once we know T 2 , we 
can obtain T\ by a Callen-Gross relationship, at least in leading logarithmic level 
(see equation |3.38|) . In the following discussion, we will leave A^ u as an unspecified 
general conserved tensor operator. Eventually we will obtain the necessary terms 
from the Callen-Gross relationship. 



3.2 The Operators 

The operators that enter into the OPE of an amplitude are composed of the fields 
involved in the interaction and derivatives. In the most general case, the deriva- 



tives can either be internal, bi-directional derivatives like the ones in equation |1.92 
or overall derivatives that act outside the fields. These overall derivatives, when 
evaluated between external states, give the difference in momentum (in momentum 
space) between the incoming and outgoing states, or, equivalently, the net momen- 
tum inflow/outflow of the local vertex of the interaction. Thus operators with overall 
derivatives are identically zero in a forward scattering process, and we do not need to 



consider them in forward OPE (see section |1.3.3| ). However, in a non- forward case as 



the one we are discussing, because the momentum flowing into the local vertex is r 
instead of zero, we have to include in the expansion new sets of operators that have 



74 



overall derivatives. 

In leading twist (the physical consequences of leading twist in this case will be 
clear later) the operators for QCD are (cf. equation |1.92|) 



= d^ n q(0)l^ + JD» n+2 ...i S MJ ?(0) (3.11) 
6™$ = d^.-.d^ F^(0)< P Mn+1 ...< 5 MJ F MJi/ (0) , (3.12) 

where it is again understood that the indices are all symmetrized. 

After taking the matrix elements between the asymmetric external states, and after 
taking a Fourier transform, the external derivatives would eventually be turned into 
factors of r M while the internal derivatives into either r M or (p+p')^ = (2p — r) M . We 
thus define the two moment variables in the non-forward 

(2p-r)-q r-q 
uj = —o , v = -zur ■ 3.13) 

Q Q 
The forward case would be the limit v = while DVCS corresponds to v = 1. 

In a QCD-Parton picture, the diagrams contributing to T Mi , are the so-called 



hand-bag diagrams as shown in figure 14. If we parametrize the momentum of the 



scattered parton as k = xp + yr (see, e.g. [16]), where x and y are two Bjorken type 
scaling variables defined by 

Qi _ i 

2p ■ q uj 

y = ^- = 1, 0.14) 

2r ■ q v 

our moment variables are related to these quantities via 
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There will be mixing among operators with the same J but different n labels. 
Under evolution in the momentum scale, which as we shall see later, is characterized 

— 2 

by Q , the internal derivatives, with which the dynamics of the process lies, can be 
turned into either overall derivatives or again internal ones, and, as stated before, 
eventually give rise to factors of either 2p — r or r. On the other hand, the overall 
derivatives, which are in essence only involved with kinematics, can be turned into 

2 

themselves only and eventually give factors of r. Therefore, evolution in Q will lead 
to mixing of these operators in only one direction in an upper-triangular fashion, 
namely, 

^ o«W»') ; o < n < ri < J . (3.16) 

However, we do still have the freedom to choose, for simplification, at a factorization 
scale /j,Q, that all internal derivatives give factors of (2p — r). Thus we can write 

(P'|6£ ( .i n ^^ (3-17) 

This is also the definition of the reduced matrix elements of OWW' 1 ). It is clear that 
from the above discussion, these reduced matrix elements will depend only on J — n, 
since the overall derivatives give (in certain sense) trivial kinematic factors of r, and 
we have 

( p '\\d^ J ^\\p) = (p'\\d {t)(J - n ' 0) \\ P ) . (3.i8) 

The choice of ( ^.17[ ) is essentially the same as taking the leading twist approx- 
imation. Analogous to the discussion in the forward case ( |1.3.3| ), the indices of the 



asymmetric matrix elements of the operators in the OPE must be made from <7 WAy or 
Pmip'm)- Aft er Fourier transformation the factors of are turned into momentum 
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factors, namely, 



And we have 



while (see |3.30|) 



WMj ^2 g2 q2 ' 



VinVlij - 2 -2 



,p-q. 



0(1). 



(3.19) 



(3.20) 



(3.21) 



Therefore the leading contribution from the operators in the OPE should all give 
factors of the external momentum, while type of contributions are small (by 
factors of =!x). It is clear that this also makes the leading twist operators dominate 
(cf. ||). Non-leading twist terms in this case are suppressed by at least a power of 
=^-. And q is now the measure of the momentum scale of the scattering. 



3.3 Wilson Coefficients and Explicit OPE 

To obtain the explicit expression of the OPE of the amplitude, we rewrite equation 
|3~T| by explicit substitution of the tensor structure operators as 

tu~)j4) VEEE^)og^(o)V--V 

Z Z J=l n=0 i=l Z Z 

+ {g m 9vp u + g^udadp - g m d v d p - g^d^) 

oo J+2 uj uj 

EEE e%{^)6%^m---—- 

J=l n=0 i=l Z Z 

(3.22) 

We evaluate the above time-ordered product between the asymmetric external proton 
states and expand it at an arbitrarily chosen factorization scale fiQ into products of 
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reduced matrix elements, corresponding Wilson coefficients, and the kinematic factors 
coming from the tensor structure operators. We obtain 

(p'\tU~)U^)\p) — 

OO J Uj 

EEE (p , ||0 W(J ' n) ||p) (Ato) i^i;i(^ 2 )(2)- J ((2p-r) • z) J -(r ■ z) n 

J=l n=0 i=l 

oo J+2 uj 

+ EEE(p , llo W(J+2 - ) lb) (M0) 

J=l n=0 8=1 

• B, u E%{z 2 ){2)- J {{2p-r) • ^) J -(r • z^C^p-r, r) . (3.23) 

Because the a and /? indices can either come from 2p—r or r when going to the reduced 
matrix elements we have included a factor of C a /3 that incorporates all the possible 
combinations in a symmetric fashion. Explicitly, it is the completely symmetrized 
sum of the following three terms: 

(2p-r) a (2p-r) p r lli r flj , 

{2p-r) a r p {2p-r) m r IXj , (3.24) 
r a r /3 (2p-r) Mi (2p-r) w , 

together with a normalization factor that we need to put in to compensate for the 
double counting of terms because, as stated before, the summation signs have already 
implied a complete symmetrization of all indices except \x and v. 

The first term of equation ( |3.23|) becomes 



OO J Uj 

E E E (p , ||0«^)|b) (Ato) (2)- J (2-C) J -"CV F%(z 2 )(p ■ z) J , (3.25) 



J=l n=0 i=l 



where we have used equation ( |2.1U| ). The second term, after a shift of labeling J+2 
J, becomes 
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OO J Uj 

EEE (p'\\0 (l)Un) \\ P ) M (g^a + g^djdp - gpa&df, - g^d,d a ) 

J=l n=0 i=l 

OO J U J 

= E E E (p'lio^ib)^) 2- (J - 2) (2 - o J - 2 c- 2 c 

J=l n=0 i=l 

• (((2p-r) At (2p-r) I/ D + g„ v {{2p-r) ■ df 

- ((2p-r)^+9 M (2p-r) I/ ) (2p-r)-d) (r-^) 2 

+ (((2p-r) M r„n + g fJ , v (2p-r)-dr-d- (2p-r) fl d v r-d - d^r u {2p-r)-d) 
+ (r /i (2p-r) I/ D + g^ u (2p-r)-dr-d 

—d fl (2p — r) u r-d — r fl d u (2p—r)-d)) (2p—r) ■ z r ■ z 

- M„ + V„) r • 9) ((2p-r)-^) 2 ) (p • z)« ^ 2jn (z 2 ) , (3.26) 

where C is the normalization factor that will compensate for the double counting. 
After using equation ( 2.10Q , it is straight forward to see that the four terms coming 



from C a/ 3(2p — r,r) all give the same value. This means that we should set C = \ in 



the above expression, which now becomes 



OO J Uj 

E E E (p'\\o {l)iJ ' n) \\p) M 2- (J - 2) (2 - c 

J=l n=0 i=l 

J-2 rp{i) 



(p.PuO + 9fiv (p ■ d) 2 - ( Pfl d v + d^ v ) p-d){p- z) J - 2 EZ, n (z 2 ) ■ (3.27) 

Therefore we have, in summary, for the asymmetric matrix elements, 
(JP'\T j,{- Z -)j v ^)\p) — > 



E EE (p'\\o^ | b) w (^H^V ^ 2 ) (p-*) j 

J=l n=0 i=l V Z 

+ 1 2 V_ 2 (PmP^ + 9AP-9? - (PA + W p-d) E%Jz 2 ) (p-z) J - 2 



79 



The amplitude T M1/ is the Fourier transform of the matrix elements with momen- 
tum q (see equation |2.18[ ). Under the Fourier transform, we have, as in the forward 
case (see section |1.3.3| , note the difference in sign convention) z^ =>■ i-g=~- There- 



t2 



fore, by taking the logarithmic derivative of q of the Fourier transformed Wilson 
coefficients, we have 



\ip^) J Jd*ze-^Eg(z 2 ) 



2^"" dtp 



- {^) J (-'?W) J ^?), (3-28) 

where ej n is the Fourier transform of Ejj n . We also have a similar equation for F. 

At the same time, by a integration by parts it is straight forward to show that, 
as in the forward case, the derivatives in the tensor structure operator will simply 
turn into factors of the q momentum, namely, => —iq^- We will drop the A^ v 
terms and concentrate on the explicit calculation of the structure function T 2 . Again 
the necessary terms generated from A^ u and thus the structure function Ti will be 
obtained from a Callen-Gross relationship once we know T 2 . We obtain the following 
results for the expression of the amplitude in the short distance limit: 

= -iJ2(p'\\d mj ' n) \\p) M (2-OT 

J,n,i 

■ (A^terms + {p^p v q 2 + g^ip ■ qf - {p^q v + q^p v ) p ■ q) 



v ■ a \ / n a 
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From equations BTO and [2.19| it is obvious 



p ■ q = p ■ q = p ■ q 



(3.30) 



Thus for the tensor structure generated by (|3.10|) in equation [3.29| we have 



Pt+Puq 2 + g^(p -qf - {p^q v + q^Pu)p ■ q 

-2 ( , (p-i) 2 p-q ( - , _ A 
= q \PtiPv + ^—=2 =r [PfiQu + Q^p^j 

C p . q Q Q 

PfiPu - -Tp-(Pnq„ ~ - P^Pu + q'^Pu + ^P^Pu) 

{p-qf / (g • g) 2 , (g • g) 2N 

^2 ' Q^Q U ^2 g2 



-2 f P-Q f , / n , (?• Q) / \ . (P- Q? ( 

q [p^Pu - -=2- + g^) + — ^- J + -t^f U» ^2 



(3.31) 



Also from equation 3.13 we have 



, ^P-Q t P-Q 
Q Q 

Therefore, we can rewrite the expression of T M „ ( p.29[ ) as 



2. . ,Z S 



T^M,q) = * / d A ze-^ (p'\T )\ p ) 



E ip'\\o mj ' n) \\p) 



J.n/, 



(Mo) uo J ^ v n ( A^terms 



P - 9 



/ (p ' g) 2 / N 

1 P^P* - ^t'kPv.Qv + + fy?" 

9 Q j 







(3.32) 



J-2 



V^rJ e« 2i Jg 2 ) ) . 



(3.33) 



Note we could not obtain the explicit expression of the A^ v terms because of its 
general form. However, there is a term of the correct tensor structure (see equation 



81 



generated from the second term. This is another indication of the Callen-Gross 



relationship that relates the two invariant amplitudes. 

Now let us define the Wilson coefficients in momentum space E as 

P ^(f^l = Hf^rj d 4 ze~^E%, (3.34) 

where we have pulled out an explicit factor of =f to make the form of E simple (see 
equation f4.60| ). We will have 



J,n,i \ 1 Q ^0 

+ ^ [pvPv - ^{P^u + q'^Pu) + {^-fq^ E%^ n (a s , -^-)J . (3.35) 

F' should be a linear combination of E and a similarly defined F. We have explicitly 
indicated their dependence on the factorization scale fio. 



By comparing with equation |2.44| we finally obtain the explicit expression of the 
invariant amplitudes, in terms of reduced matrix elements and their corresponding 
Wilson coefficients, as the following: 

HT, = - £ (p'llOW^Ib)^^^^'?^^, 3) (3.36) 

J,n,i A*0 



TH^CT2 = ~ E (^iiO (i)(J,n) lb)( M0 )^ J ^ ri ^2 I n(« S > ^) • (3-37) 

J,n,i "0 

We will concentrate on the discussion of T 2 from now on. We always regard 
T\ as being obtained from T 2 from a Callen-Gross relationship that can readily be 
extracted from the above discussion: 

Ti = ^|T 2 . (3.38) 
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Chapter 4 

Renormalization Group Analysis 

In this chapter, We perform a renormalization group (RG) analysis of the non-forward 
amplitude. We write down the renormalization group equation (RGE) of the oper- 
ators and their corresponding Wilson coefficients of the operator product expansion 
(OPE) of the amplitude, and find the formal solution to the RGE. Then we compute 
the evolution kernels that will generate the anomalous dimensions and discuss their 
properties. Finally we will calculate explicitly the lowest order Wilson coefficients. 

4.1 Renormalization Group Equation 

As stated in the last chapter, we will concentrate on the invariant amplitude T2, and 
Ti can be obtained from the Callen-Gross relationship ( p.38| ). T 2 is analogous to the 
structure function F 2 (1/W2) (see section |1.2.2| ), however, the analogy is not complete 
because of the subtleties in the discussion of dispersion relationship later on. 
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Recall that T2 is given, after OPE, as 



P-Q 

M 2 



1 - ct 2 = - e mo^ n) \\v) { ,y^ n Ef^ n 

J,n,i 



S ' Mo' 



(4.1) 



It can be regarded as a double distribution function. In the usual language it is the 
series sum of "double" moments in terms of the two new moment variable u and v 



that we defined earlier (see equation |3.13|) : 



T 2 = Y,uJ J - n v n T { 2 J ' n \ 



(4.2) 



J.n 



with the double moments T^ J ' n ^ defined as 

:P-9rp(J,n) 



1 - C— L T ( 

S M2 2,(Q,/k>) 



i Mo 



(4.3) 



We will eventually analytically continue in J, but leave n as discrete. The dis- 
cussion of analyticity properties of T2 is almost identical to that presented for the 
forward case (see section |1 . 3 . 3 ) with minor adjustments. Explicitly writing out the 
spin sum of T Mi , we have (see equations |2.16| and |2.18|) 



T^(p, q',q) 



*ze-^{p',s\TU-^j u ^)\p,s) 



2 

= iJ#ze^*{jt\TM-z)j v {Q)b) 

= i W d'ze^{e(z )(p'\j fl (z)\r}(r\j l/ (0)\p) 

r J 

+9(-z )(p'\j v (0)\r)(r\j,(z)\p)} 
= i W d 4 ze^ z {e{z Q )e^'-^ z {p'\ 3 M\r){r\ 3u mp) 



+e(-z )e 



-i{p- Pr )-ZI I 



- z (p'\jM\r)(r\jM\p)} 



(4.4) 
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where we have used the translation operators (see equation 1.18 ) in the last step. 
Performing the spacial integration yields a 3-d 5-function and we have, similar to the 
forward case, 



T = % 



(2n) 3 J2 I ^o(5 3 (g + ^ -Pr)e 4(90+p ' ^ o) - z ^(^o)(p , |^(0)|r)(r|j,(0)b) 

+(5 3 (g-p + p r )e^- po+ ^>^^(-^)(j9l^(0)|r)(r|j M (0)b)) 
_ (27r)3 y f S 3 (q + P'-Pr)(p'\j,(0)\r)(r\3v(0)\p) 

r\ QO+P'o- Prfl + it 

S 3 (q-p + p r ){p'\UO)\r)(r\j„(0)\p)\ 



(4-5) 

Qo ~ Po + Pr,o -ie J 

Again by looking at the poles in qo we can discuss the analytic properties of the 
amplitude. The first term gives rise to poles at 



QO + Po = Prfl = E r 



(4.6) 



which means that (cf |1.3.3|) 



?'\2 



(g + PoY = K = M 2 +f r = M 2 + (q + p') 



(4.7) 



where in the last step we have used the spacial ^-function. M r is again the invariant 



mass of the intermediate state r. Thus we have, recall equations 2.1 and 1.1 



M r 2 = {q + p'Y = {q + pY = q' 2 + 2p-q 



Rewrite q' 2 as 



/2 1 / 12 , 2\ 1 / 2 I2\ 

Q =2(9 +9 ) - -? ) 



and use equations 2.S and 3.13, we arrive at 



(4.9) 



Ml 



(q' +p) =Q +{2p-r)-q 



q 2 {l~u) 



(4.10) 
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and thus the poles in to are at 

M 2 

LU r = 1 + — f . (4.11) 

Q 

The second term of the spectral expansion of the amplitude gives poles at 

9o ~Po = -Pr,o = ~ E r (4.12) 

and after similar discussion results in 

M 2 = (q — p) 2 = q 2 — 2p ■ q = q 2 — (2p — r) ■ q 

= t(l+u), (4.13) 

which means that there are also poles in uj at 

M 2 

0J r = -1 - — r ■ (4.14) 

Q 

We see again that the momentum scale is set naturally by Q = —q 2 . 

It is clear that in this non-forward case, also has an analytic circle of unit 
radius in the a;-plane and has a Taylor expansion in u. The discussion on the dis- 
persion relationship for T 2 then follows almost identically to that of the forward case 
(see section pL .3 .3| ) 

We can indeed use the dispersion relationship to analytically extend the discussion 
from the region of uo < 1 where T^ is analytic but kinematically unphysical to the 
physical region of the Bjorken limit where < (cu) _1 < 1 in the same manner as in 
the forward case. And we find that the inverse of the expansion W?2 



is 



T^ } = /-/-(V-(VT 2 (Q 2 ,^,) 

Jc UJ Jc h> UJ V 

r°° du f dv , 1 , 



Jl UJ Jc is UJ v 
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where c is any contour in the 00(1/) plane that encloses the origin and W2 = ImT2. 
Note that in obtaining the second equality we have used the fact that T 2 is even in 
uj. Also note that the moment index of uj in this case is, in our notation, J — n. 

T2 now as a double distribution function is given by a Mellin transformation: 

T 2 (Q 2 ,uj,u) = / !iVe (J - n)log %' 1 T? ,l) (4.16) 
Jc 2m n=0 

with the contour c lying parallel to the imaginary axis in the J — n plane and to the 
right of all singularities in that plane. 

Under renormalization, the operators scale according to a renormalization group 
equation 

^0^>=EE^^ (4.17) 

" n'=n i' 

where 7 is the anomalous dimension matrix which acts in the product space J <8> uj. 



Because the mixing of these operators is such that n — > n' >n (see equation |3.16| the 
matrix 7 is upper triangular in the J dimensional n-space. 

Since T 2 is a physical quantity it should not depend on the factorization scale. The 
Wilson coefficients obey a similar renormalization group equation thus cancelling out 
the scale dependence of the operators: 

^ 2 4r,E%{a s {^ ?) = - ± til (4-18) 

ufJ> H> n / =0 v 

We can write the solution to this RG equation as 

E%{* s {v),% = E E^3(« S (Q),1) Mtf n (a(Q) 5), (4.19) 

f 1 n'=0 V t 1 

where M is a path ordered exponential of the anomalous dimension matrix, formally 
given as 



<Ua S (Q),®) = (veM- /J7(«,(A 2 ))^) ) . (4.20) 
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Therefore by substitution the double moments are given by 

= -± E« 2 ,,K(Q),1) M^ n (a s (Q)A) (P\\d^\\p) (p0) . 

n'=0 i'i r 1 ® 

(4.21) 

In principle, for the high energy scattering process, if we know the reduced matrix 
elements at a given momentum scale //, we can calculate the Wilson coefficients 
order by order in perturbation theory at some high momentum scale Q, at which the 
scattering actually takes place and perturbative QCD is valid, then use the anomalous 
dimensions, also computed to some fixed order in perturbation theory, to evolve along 
the momentum scale to \i so that we can use these reduced matrix elements and obtain 
(theoretical predictions of) the invariant amplitude, the amplitude and eventually the 
cross section. 

The matrix elements can come from two sources in general. They can either come 
from a first-principle calculation in lattice QCD, or from a high energy experiment. 
The former is phenomenologically not quite feasible yet. For the latter, we actually 
measure the cross section of the high energy scattering at a certain high momentum 
scale, use the above mentioned procedure in reverse and phenomenological models 
to extract the matrix elements, and then by RG analysis we can predict the cross 
section at any other high momentum scale where factorization is valid and compare 
the results with actual experimental data at the scale. 

In the next two sections, we calculate the anomalous dimension (or rather, evolu- 
tion kernels from which they can be extracted) and Wilson coefficients in the lowest 
non-trivial order. 
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4.2 Anomalous Dimensions 

We will use the light cone (LC) gauge and compute the lowest order anomalous dimen- 
sions in terms of light cone variables. In the LC gauge, the lowest order anomalous 
dimensions of the quark and gluon operators are generated by the usual triangle dia- 
grams. These diagrams are shown collectively in figure |15[ The graphs where we have 
more than two quark/gluon lines meeting at the vertex are zero because of the gauge 
choice. The price we pay is the presence of extra end point colinear divergences in 
the integration of light cone variables, which, as we will see, will cancel out eventually 
after we include the self-energy graphs. 

We define the momentum fraction variables (which are the moment variables) as 

W = ^, 1/ = ^. (4.22) 
k 1+ k 1+ 

The conventions (cf. ref. ]T4]] ) we use for the quark and gluon vertices are 
{ q J ' n)i+ \k,k-r) = 7+ (2£;-r)^- 1 r™ 

Of n)aP {k, k-r) = 2g aP n- kn ■ (k-r)(2k- r)f^r£ , (4.23) 



where n is the LC null vector (see equation |1.119|). 



4.2.1 Quark-Quark Anomalous Dimension 



The detailed version of the first diagram in figure 15 is shown in figure 16, where a 
and P are Lorentz indices; a, a', b and b' are color indices for the quarks while i labels 
the gluon line. In contrast to figure |9] of the forward case, the gluon is not on-shell. 
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Using standard Feynman rules and the LC gauge we can write down its value as 



i.b 



7+' 1 = (fc) 2 £(lML0 / 7^^ 1 ^0^ + \k,k-r)- 1 p 



iD a ^(k 1 -k) 



(2tt) 4 ' a jt- / 



(h - kf 



(4.24) 



which, after the color algebra (see section |1.4.3|) gives (where we have taken the 
incoming and outgoing quark to have the same color index) 



d 4 k 



lt q l = (i9?C F ^-^-J—O^ik, k - r)^p 



iD a < 3 (ki - k) 



(2vr) 4 / 



$ " (h - k) 2 



(4.25) 



where the quark vertex O q is given in ( 4.23|) and the LC gluon projector (numerator 



of the LC gluon propagator) (see equation [1.121 



D a/3 (k) = g a/3 - 



p _ n a kp + k a np 



n ■ k 



(4.26) 



Again n is the LC null vector. Writing the integral of the loop momentum in terms 
of light-cone variables (see section |1.4.3|) we get 



.+■1 _ : !l ' r ' I .-a 



d k dk+dk- 



(2k ~ry- n ^rlA(k,r) 



(k 2 + ie)((k - r) 2 + i€)((h - k) 2 + ie) 



(4.27) 



where 



A(k,r) = 7a (#- /)7+ tflpD^ih - k) . 



(4.28) 



We perform the integration first, in terms of the poles of the integrand, which 
come from the denominator factors 



k 2 + ie = 2k + k- - k 2 + ie , 
(k-r) 2 + ie = 2(k - r)+£;_ - k 2 + ie , 
(h-k f + ie = 2(k - h) + k_ - k 2 + ie . 



(4.29) 
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However, the positions of poles in fc_ are different when k + is in different regions. 
More specifically, there are four regions we need to consider: 

f ) k + < 0, by setting the denominators to zero we find that the poles are at, respec- 
tively, 

k - = w + +ie ' k ~ = w^ +ie ^ k - = 2(k-k l)+ +i€ - (4 - 30) 

All three poles lie above the real axis of the complex plane (see figure pj] (a)). We 
can complete the contour of integration of k_ in the lower-half-plane and clockwise, 
enclosing no poles and thus the integration gives zero value. 

2) < k + < r + (< ki + ), the poles are at, respectively, 

k - = W + ~ ie > k - = W^r)- + +l ^ k - = 2{k-k l)+ +ie - (4 ' 31) 

The first pole is below the real axis while the other two are above it (see figure |TT| (b)). 
We can complete the integration contour again in the lower-half-plane and clockwise, 
picking up the pole at k- = — ie = k^ . The value of the integration would then 
be the residue of the pole with a factor of — 2iri. 

3) r + < k + < ki + , the poles are at, respectively, 

t 2 k 2 k 2 

k^ = -= ie , A;_ = ~ — r ie , k_ = - ~ - . h ie . (4.32) 

2k + ' 2{k-r)+ 2(k - k 1 ) + v ' 

The first two poles are below the real axis while the third one is above it (see figure 



17\ (c)). We can complete the integration contour now in the upper-half-plane and 



counterclockwise, picking up the pole at k_ = 2 (k^-k 1+ ) + ^ e = ^ • The value of the 

(2) 

integration would then be the residue of the k_ pole with a factor of 2ixi. 
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4) k + > ki + , the poles are at, respectively, 



JL 

2k A 



— ie , k- 



k 



2(k-r) + 



— ie , k. 



k z 



2(k - ki). 



— ie . 



(4.33) 



All three poles lie below the real axis of the plane (see figure |I7] (d)). We can com- 
plete the contour of integration of k_ in the upper-half-plane and counterclockwise, 
again enclosing no poles and the integration gives zero value. 

Thus the k + integration should be performed in two regions: 

,j2Cf ff k ( r w+ r*. ^ ' ,ik - < 2fc - r ^ r * r » 



7 



gg 



—i 



(2vr) 4 

.CY s Cf 

dk. 



o 



r+ 



(k 2 + ie)((k - r) 2 + ie)((Jfei - k) 2 + ie) 



dk 2 



dk + + dk + \ (2k - r 

o Jr+ J 



J^n—l^n 



A(k,r) 



(2k+k. - k l + ie)(2(k-r) + k_ - k 2 + ie)(2(k-k 1 ) + k. - k? + ie) 



(4.34) 



Evaluating the contour integral of fc_ at the poles specified above, we have 



7 



+,i 
gg 



-—i 



(2tt)2 y dtL 

f k l+ 
+ I dk + 



dk A 



-2m) 



(2k-r) J - n - 1 rlA(k,r) 



o 

(2?rz) 



2£;+ \ k (2(fc-r) + jfe_- J fc 2 )(2(A;-A; 1 ) + A;_- A; 2 
/ (2k-r) J ^- 1 rlA(k,r) 



k-=k 



(i) 



r + " t "2(A; — V(2fc + fc_-r)(2(fc-r) + fc_- £ 



(2) 



(4.35) 



Substituting the pole values into the denominators and simplify them first, we get, 
for the pole k^ = ^- , 



(2(fc-r)+fc_- fc 2 )(2(Jfe-Jfei) + A;_- A; 2 



fc_=fc 



(i) 



^ 2 (fe-r) + _ fc2 \/ 2 (fc-fci) + _ fc2 \ , , ,,, 



mi _^zik)( 1 + ^-k). 



k A 



k A 



1 fe. 



k A (k+ - (k-r) + )(k+ + (k-r)+) k 4 r + k 1+ ' 



(4.36) 
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while for the pole k 



(2) 



2(fc+-ki+) ' 
2 >,.),/. L 7„2\ 



(2k + k_- k 2 )(2(k-r) + k_ - k 2 



2 ^+ (2\//2 r ) + i 2\ /;2\2 / i , ^+ 1 / i , — r )- 



1 (*i-*0+ 1 {h-k) + 



k ((h-k) + + k+)((k 1 -k)+ + (k-r)+) k A (k-r)+k 1+ 
Thus we arrive at 



(4.37) 



^ — r~TT r- dk +( k i ~ k)+A {2) (2k - r) J r^ l rl J , (4.38) 



-2fci + (fci — r 

where A* 1 ) and are A(A;,r) evaluated at the poles k^ and k^\ respectively. 

We are interested in the logarithmic divergence in the graphs so in computing the 
Dirac spinors we keep only the leading powers of k. This leading power turns out to 
be quadratic, which means we will keep in the following calculation of and A^ 
only terms that are proportional to fc 2 . 

Defining u> and v as before in (|4.22|) , we have, at the poles, 



}^ = ^- = -=— & k 2 = 2k + k^-k 2 = 0, 
2k + 2k 1+ uj 

k 2 — k 2 k 2 

k (2) = « = K & p = « 

2(A; + -Jfei+) 2Jfei+(l-w) 
Recall that A(fc, r) can be written into two terms: 

A{k,r) = la (}l- /) 7+ h:aD aP (h-k) 

= la HpD aP {kr -k)-j a h+ ^D a(3 (h - k) 

= At -An (4.39) 
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Using the anti-commutation relationships of the Dirac matrices in light cone 
variables (see equations |1.115| and [L.117| ) we can reduce the first term (Aj) to 



la ftlpD^ih - k) = 2k +la ^^(h - k) - k 2 lal+1 pD a \k x - A01.4O) 

At the pole k_ , because k 2 = 0, we have 

At = 2k +la $ lp D^(k x -k) = A { p . (4.41) 

Expanding D al3 (ki — k) and again using anti-commutation rules and the fact that 
k 2 = we have 

A? = 2k +la ^ a - - ^ + ( 7+ flfr- + (ft- # ^ 7+ ) 

fci + (l - 

= -Ak + j{--^L( 1+ k' ft+ ft - 2 1+ k 2 ) 

= -4k + (k +1 „ + k„ 1+ - k ■ 7 ) - t^(7+ ¥ h+ fa ^7+) • (4-42) 

— l — UJ 



As stated before, we want to keep only the leading power in k. It is clear that in ( |4.42|) 



the leading power of k is quadratic and can only come from a term that contains the 
factor k- . In addition, since we are computing the quark-quark anomalous dimension, 
only terms that have similar vertex structure as the original quark vertex O q in ( 4.23 ) 



contribute. Thus, we keep in ( 4.42| ) only terms that are proportional to 7+A;_. By 



expanding the dot products and inspecting the result, it is straight forward to see 
that only the second term in the first parenthesis survives and we have 

A ( p = -Ak+k_ 1+ = -2k 2 <y + , (4.43) 



where we have also used ( |1.116| ). 
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The second term of (|4.39p can be similarly reduced by keeping only the terms 
proportional to r y+k_: 

4V = la /7+ ftli3D al3 (h - k) 

= 7 a h + ^ - 77 * (7+ /7+ Hh~ + (h~ t) /7+ ^7+) 
[ki - k)+ 

= - tu 1 (2r +7+ (* ft - A: 2 ) + (ft /- jf /)2fc+ 7+ ) 
(«i - «)+ 

= ~ (fc 1 -fc) + (2r+7+(7 - fc+ ^ " °) + (°- # /) 2A; +7+) 

-(0 - 2A; + # / 7+ ) = - — — — 2A; + fc_r + 7 + 7_7 + 



(*!-*)+ + ' ' (fci - 

2/ " i 2 7+ = T^-fe 2 7 + • (4-44) 



(&4 — fe) + 1 — OJ 

Note some of the zeros above do not mean numerically zero but rather the contribu- 
tions of those terms are zero. We use zero (0) this way through out the discussion 
unless otherwise stated. 

Combining the results we have 

= -2A; 2 7+ (1 + — ^— ) . (4.45) 
1 — OJ 

The evaluation of goes parallel but is a little bit more involved since we 
no longer have k 2 = but rather k 2 = — ■ We outline the main steps in the 
following: 

= 2k +la Jt^D^ih - k) - k^j+jpD^ih - k) 
2k 

= 2k +la ^ - + ( 7+ Jf(fr- t) + (h~ JQ 

fcl+(l - OJ) 

-k 2 lal+1 a + —. ^ r(7 + 7 + (^i- ¥) + (h~ 07+7+) 

fei + (l - OJ) 

= -4k + ft - ft h+ h - ^1+k 2 ) + 2k 2 j+ + . 

1 — OJ 
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Using identical arguments as those in calculation of plus the fact that now the 
leading contribution (of A; 2 ) can also come from a k 2 factor in addition to we 
arrive at 

Af = -4k + k. 1+ --^{-2 1+ k 2 ) + 2k 2 1+ 

-k 2 , 4 ^ „x , k 2 , 
= -4fc + 0/ , ~, , 7+ + (^ + 2) 7+ ( 



A; 2 



„2. 



2(ib + - ibi + ) '" r v 1-cj y (1-w)' 
7+(2cj(1 - w) - 4a; - 2(1 - w)) 



2 r7+ „ , . 2 



(l-^) 2 

On the other hand, 



(1 + cj 2 ). (4.46) 



4? = 7a /7+ ^lpD a/3 (k 1 - fc) 

= 7. h + - (fcl ^ fc)+ (7+ /7+ ft ft- + (^i- /7+ ^7+) 

= - tv 1 i,\ ( 2r +7+(^ ft " A: 2 ) + (ft /- ? /)2A; +7+ ) 
(«i - fc)+ 

= ~ {k 1 -k) + {2r+1+{l - k+ h ~~ e) + ( °~ V /)2A;+7+) 

= z k) + ^ 2r+1+k2 ~ 2A: + /c - r + 7 + 7 - 7 +) 

2r, -k 2 Ak+ k 2 



"7+7^ \ + T, T^7+^+ 



(fci-^+'^a-w) (^1-^)+ " + 2(fc+-fc 1+ ) 

- 7+ U + (l-«) 2 + ^ + (l-o;)V 

2 ^ 7+ (!+")■ (4.47) 



Again combining the results we have 



A^ = -^y 2 (l + . 2 -u(l + .)). (4.48) 
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Thus we have 
7+' 1 = ~C F 



f ^ 7+ Jfe£; V (- r du (2u - u) M uj (1 + - 
J k \u Jo 1 



+J- fdu^-v)^ I±^0±^A . (4.49) 

1—UJu 1 — LU J 

The a; integration is divergent at u; = 1. This divergence is cancelled by the 



self-energy graphs, the value of which in LC gauge can be readily taken from [47] as 



Z F (x) = 1 + ^ V2/0 - 2 log |x| + ?) , (4.50) 

Z7T e z 



where Jo is the colinear divergence 

- 1 dz r 1 dz 



'0 



(4.51) 



10 z Jo 1 — z 

This divergence can be viewed as essentially an artifact coming from the choice of 



using light cone calculations. For a detailed treatment please refer to |47|. Here we 
simply state and use their conclusions. 

Zp depends on the longitudinal momentum fraction x. So in our case the self-energy 
contribution from the k\—r line should have x = n '^ 1 fc ~ r ^ = 1 — v while the k\ line has 



simply x = 1. Adding |(Zp(l) + Zp(l — v)) to (449) and identifying the logarithmic 



divergence / = ~ we will obtain the expression of ^ qq as in (|4.53| ), namely, 

7# = ^C Fl+ k(-\ n f du{2uj-v) M (®{v-uj)-uj{l + -^) 
2tx Jo \ v 1 — io 

+ 9 ( M _ y) _L '+"'-Ki + ^) + , (1 _ w)( _ 2/o _ log(1 _,) + »,) . 

1 — v l — ^ / 

(4.52) 

Note there are no other diagrams contributing in LC gauge. Recall this is because 
the gauge choice makes diagrams that have more than two lines meeting at the vertex 
vanish. 
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4.2.2 Anomalous Dimension Kernels 

While the details of the computation of the other three anomalous dimensions are 
shown in the Appendices, we list the result in the following: 



Iqq 



- ^C Fl+ k(- l u n C duo(2uo-u) J - n - 1 (e(y-u)-u}(l + -^—) 

Z7T JO V V 1 — iO 

+ q( u - v )^- 1+u \- u{1+u) + s(i-u)(-2i - mi-,) + 5) 

1—vl — u 2 

7$ = ^C Fl+ k(- l u n f l du(2uj-u) M (q^-uj)-^ 2 ^) 
a 2n Jo \ v 



iqg 



•39 



— -2g lj ki,v n j 1 duo(2uj-v) J ^ l (-Q(v-uj)-uj(2uj-l-u) (4.53) 
27T 2 Jo V V 

+Q(uj-v)-^(u 2 + (1-u) 2 -uji>)\ , 
l-u J 

?^2g ij k(,is n FdLO&u -p) m (q(u-uj)-(- (Alu 3 - lu 2 + 4u) 
2tx Jo \ 2 v 

^+uj 2 -2uj 2 u 2 Q(u-v) , 2(l-u+u 2 ) 2 z/(l+u; 2 )-2(l+a; 3 ) N 
+ ; 3(w +w))) + — 

1 — C<J 1—1/ 1 — U 1 — LO 

+ (l-^)5(l-a;)(-2Jo-log(l-^) + ^-)J . 

We have included virtual corrections coming from the self-energy graphs that will 
cancel the colinear singularities at end point in the u integration (J terms, see equa- 



tion [4.51| for the definition of 7q) an d give correct constant terms in the anomalous 



dimensions. is the usual step function and 

b = 11 - %i f (4.54) 

is the leading coefficient of the QCD /3-function. 
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4.2.3 Comments on Anomalous Dimensions 

Some comments are in order. It is straight forward to show that in the forward 
limit where v — 0, equation [4.53| reduces to the conventional Altarelli-Parisi splitting 



functions (T^|. However, because of the non-forwardness, there is no simple probability 
interpretation for these evolution kernels as splitting functions. 

The reason we have two terms for each triangle graph is that for different k + 
integration regions pinching of the pole is different (see section |4.2.1| ). Thus there 



is no straight forward optical-theorem type dispersion relationship between the cross 
section and the imaginary part of the amplitude. Nonetheless, we can still do analytic 
continuation in J and relate the matrix elements to "double parton distributions" ]12| 



which now may not have a direct probability interpretation. The second terms of our 
7 differ from those of [^| only because the factor n- kn- (k — r) in our convention for 
the gluon vertex (see equation ^4.23| ) is not included in the definition of their gluon 
vertex. Once we take this into account and use the same vertices, their results are 
identical to the second terms of the corresponding kernels of ours. [] 

After we perform the uo ( i.e. k + ) integral, we can in principle obtain the cor- 
responding anomalous dimension matrix in the moment space. At first sight, the 
(1 — u)^ 1 factor in the second terms may generate a series of infinite sums over pow- 
ers of v, which will spoil the locality of the vertex thus invalidate the operator product 



1 Simultaneous to our work, Blumlein et al |4§| also calculated the evolution kernels of what 
they call "twist 2 light-ray operators for unpolarized and polarized DIS" . It would be useful to 
make a detailed comparison to see whether the seeming deviation between our results is a nontrivial 
disagreement or mainly due to different notations and approaches. Because of these differences in 
notation and approaches, together with the fact that their main results egs.(15) — (18) were stated 
but not derived, we are not able to make such a comparison at the present time. 
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expansion. Detailed calculations show that because the lower bound of integration is 
now v instead of zero, there will also always be at least one power of (1— v) coming out 
from the integral. Furthermore, the highest powers of v cancel completely between 
the first two terms of each 7, leaving the highest surviving v powers to be v J ~ x for the 
quark sector and v J for the gluonic sector, as they must be to make the OPE valid. 
All the Iq and log(l — v) terms that might potentially spoil the operator product 
expansion also cancel completely between each triangle graph and its corresponding 
self-energy graphs. 

Because of the mixing among different n moments, as well as between quark and 
gluon sectors, it is very difficult to read off the anomalous dimensions of the mixing 
between two operators (with same J) of definite n moments. Q However, the form 
of the anomalous dimension matrix simplifies greatly when the high energy limit is 
taken, and we will be able to make connection between this general formalism and 
the conventional leading logarithmic approximation (LLA) analysis results. 



4.3 Lowest Order Wilson Coefficients 

The lowest order Wilson coefficients can be calculated in perturbation theory from 



tree level Born diagrams as shown in figure [18], where we have explicitly included the 
cross diagram. 

Using standard Feynman rules we can write the value of the Born diagrams, after 



2 After the completion of this work, through private communications with Professors X. Ji and A. 
V. Radyushkin, the author learned that by adopting a basis of linear combinations of the operators 



in 3.11 using Gcgcnbauer polynomials as coefficients we might be able to explicitly diagonalize the 



anomalous dimension matrix in the moment space. 
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spin average, as 



T (o) 



1 Y^r ^77 / a ( A') 

oLK U s{P) Tm 7 T a2, • 7* 



+7, 



(p — g) 2 + ze 



7m )U s (p). 



(4.55) 



By using equations |2.10| , |4.10| and |4.13| we find 

1 



T (o) 



where 



-let 



^ 2q 2 
1 



E^(p)7m(^+ ^)7^.(P) 



1+w ... 



v 1 - w 

E^(P)7,(i>- ^)7^ S (P) 



-ze 



oo 

X><(Tr(i> 7/ ^+ + (-l) l Tr(Mft- fiyj) 



q 



1=0 
oo 



E w '(pm(p + 9% + pAp + q')n - 9v»v ■ (p + q')) 



1=0 



+(- 1 Y(p f i(p - q) u + Pv{p - qOai - sw? • (p - <?)) 

oo 

■ E u l ((p^Pv + PvPv + Vydv + Vvin - 9»p ■ q) 



-9 g2 



i=0 



+(-i)'(p^ + - PpQ» - v v % + 9»»p ■ q)) 

oo 

E u; 2l (2{p IJ ,p u + MV) - p M rv - r^p u 



'9 g2 



Z=0 



+w(Pn(tf + q)u + (q' + q)nPv ~ 2g ftl/ p ■ q)) 



9 Jf 

W, 1=0 



(4.56) 



CV = ( 2 P - r )» r » + r n( 2 P - r )v + Mp^ + V* - ^P • ?) • 



(4.57) 



By similar derivations that lead to equation ( |3.31| ), as well as the definition of the 
moment variables (see equation |3.13| ) we find 
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a 



2(2 




2(2 


- QipuPv - 


2(2 


-<){,*, 


2u 


[p-q 



p-q , 
t 1 

p-q 



{p ■ q)'' 



(p^qv + q'uPu) + (—^-fq'uqu + 



(p-q) 



_ 2 \f^v < Ififuj 1 V - 2 / '/ 



,p-q,2 , 

q F 



— (a 

2 q2 , 



Thus we obtain the value of the Born diagrams as 



(4.5* 



T (o) 



i=o \ q q 

1 ( P'q? i / \ i (P - q \ 2 1 
+ — I PyPv - -=2-{p^qu + q^Pv) + (-^-) q^ 



QuQu, 



p-q 

Comparing to (|3.35|) we obtain 



q 



q 



(4.59) 



Aie q , J odd, J > 3 , 



J, n otherwise . 



(4.60) 



Thus for each J value, only the Wilson coefficient of n = has a non-zero value 
at the tree level. This means that at leading order there is no dependence on the 
non-forwardness in the Wilson coefficients. 
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Chapter 5 



High Energy Limit 



In this chapter we will solve the renormalization group (RG) equation for the Wilson 
coefficient (|4.18 ) in the high energy limit, and discuss the relationship between our 



general results and that of the usual forward double leading logarithmic approxima- 
tion (DLL A) 0. 



5.1 Reduction of the Path Ordered Exponential 



Since we have computed the lowest order Wilson coefficients (|4.60|) , and the reduced 
matrix elements are regarded as input from experimental data, we need now to evalu- 
ate the path ordered exponential M^ n (a(Q) , — ) in equation (^.19|) in order to obtain 



the double moments T{- J !^ . and then the invariant amplitude T 2 . 

2,(Q,Ato) v 1 



Recall that formally we have ( [4.20 ) 



M&(«.@),-) = (pexp( 



7(«,(A 2 )) 



d\ 2 



(5.1) 



n'n 
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thus to explicitly solve the RG equation we need to diagonalize the anomalous dimen- 
sion matrix in both flavor (i) and moment (n) indices. However, in the high energy 
limit the situation simplifies. 

Similar to the forward case, after analytical continuation in J, at high energy, 
the dominant contributions to T2 come from the leading (right most) poles in the J- 
plane, which, as we will see, means that the diagonal elements in the moment space 
(n space) of the anomalous dimension matrix dominate the evolution. And because 
the gluon anomalous dimension has the right most pole (one unit higher than that of 
the quark sector) (see equation |5.4| ), together with the fact that gluons have bigger 
color charge, in flavor space (J, space) at high energy the gluon anomalous dimension 
dominates. 

However, gluons only have color charge and thus can not interact directly with the 
color-neutral virtual photon. The scattering must happen via a quark loop, which 
means we have to force at least one gluon-quark transition at the end of the evolution. 
Thus in each term of the expansion of the path ordered exponential (|4.20 ), all the 



factors of the product of anomalous dimensions are 7" except the first one, which 
should be 7^ due to the quark-gluon transition. 

Therefore we need only to evaluate M between a quark and a gluon state and thus 
collapse the two flavor sums in ( |4.21| ) to i! = q and i = g. Formally, if we now label 



j u> as the anomalous dimension matrix in the moment space, by expanding the path 
ordered exponential, we have 

_ Q il=q>i=g ( f Q 2 ^ f ^dX\\ qg 



Mj,« B (a.(Q),^) l ' =J ± =9 (VeM-J 2 7WA 2 )) AL ,,, 

" V ^ / n'n 

(00 1 „q 2 1 J\2 

1=0 L[ -V k=i At 



n'n 
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e L ^f(-iM\m /; ^(-7K(a 2 h ))) • .. 

2 ^(-7(a.(A?)))j 

1 / n'n 

(oo .q2 j\2 „A 2 W\ 2 

E / ^(-7 W MA?))) /; ^M"Mxh))) ■ 
i=o J v A i J » \-\ 

d\ 2 \ 



\i 



^0 A ■'^0 A 



- 7 "(«,(A 2 ))Pexp(- / -^7 59 (« S (A /2 ))) 



(5.2) 



which can be rewritten as 



M^(a s (Q),^)= -/ -w7 9S K(A 2 )) ^exp(-/ — 7 99 MA' 2 ))) . 

/^O V ^0 A' / n /, V J & X Jin 

(5.3) 

The leading pole terms of the relevant anomalous dimensions can be parametrized as 
(suppressing the J label), after integration of (|4.53|) , 



r? 1 4- r" 

7^(A 2 ) = -^a a (A 2 )^,, with >C = j n ' ,U i ( 5 - 4 ) 

7^(A 2 ) = -^« S (A 2 )7^, with 7^ = J _ n ":" 1 • ( 5 - 5 ) 
where < n < n' < J for 7 nn , and more importantly 

<£U„ = • (5.6) 

This means that only the diagonal elements have leading pole contributions. 

Because of the upper triangular structure and the leading pole positions (recall that 
the right most pole dominates) of these 7 matrices, we have, for m < I < n, 

y / 1ml lln Imn Inn = Iran flrl ' (5-7) 
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where we have labeled the diagonal terms of the anomalous dimensions (leading pole 
terms) a n . This means in a product of 7 matrices we can, in leading logarithmic 
approximation (LLA), take only diagonal entries in all the factors except for the first 
one (see equation |5.11| ). 



On the other hand, because we have the QCD running coupling (see, eg 

47T 

6 log"(A 2 /A 



where A = Aqcd is the QCD scale and 60 the leading QCD /^-function coefficient (see 
( [4. 54] )), the integration of the anomalous dimension at high energy gives 

±C A _ QQ ^rfi og (A'VA 2 ) 



-7 



7T ' J,l log(A' 2 /A 2 ) 

= A(A 2 )7 OT , (5.9) 
where we define the usual double-log factor 

Thus the path ordered exponential of r y" (the second factor on the right hand 



side of |5.3| ) becomes 



= (l + ^7 + ^ 2 77 + -)in 
A 2 

= (l + ^7/n + ^f a n7/n + --) 

= l + (A + ^a n + ...)^ 



I+ J_ (e ^K_ 1)Tf , 



a?, 



}_ e A(\ 2 )ai ^99 + non — leading terms . (5.11) 

On 
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Therefore, (15. 31) reduces to 



\ J Ho A / n'l 0,71 





rQ 2 dX 2 


(- 


Jfi A 2 




rQ 2 dX 2 




U A 2 


a± 


rQ 2 d\ 2 


2w 


U A 2 



a s {X 2 ) e A ^^ m ^-j 



1 

nJX 2 )e A ^^, n , (5.12) 



where in the last step we have again used (|5.7f) . 

The n! summation in the expression for the double moments ( |4.21| ) also collapses 
in LLA to the diagonal element of 7 g9 , and we have the final expression of the double 



moments as 



n'=0 i'i 

= - E Ef_ 2 ^a s (Q)A) KIMQ),-) (P\\0 {9)ij ' n) \\p) ( , o) 

n'=0 

= - t £&,„.K(Q), 1) £ /f ^<MA 2 ) e^'*" 

n'=0 27F J ''o A ~ 

= -^(« s (g),l)^/J^(A 2 ) 

■ri(p'\\o^\\p) { , 0) 

= -E { ; q ! 2 , n (as(Q),l) ¥WO^ J ' n) \\p) M 

n f [Q 2 dX 2 2 a(x*) 1 

' o / , -TT«- A 2 eJ -n-i_ . 5.13 

27T j/ig J — n 

Now the invariant amplitude T 2 ( J4.16 ), after analytically continuation in J, be- 
comes 

- PyWT 2( <f , «, *) = / £ J>" /J f 1) 

. (p '||oW«»)||p> w J_e X p((7- n )loga, + - 7 ^ T ) . (5.14) 
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In the high energy limit we can make a saddle point approximation for the J 
integration because the structure of the exponent in the expression. Define 

A 



f(J) = (J-n) logcu + 



J-n-V 



(5.15) 



we obtain the first and second order derivative of / as 

A 



f(J) = logc 

/V) 



2A 



(J-n-iy 



(J-n-lf 



(5.16) 



Setting f'(J) = we obtain the saddle point J s as 



J s = 1 + n + 



logo; ' 



(5.17) 



and at the saddle point we have 

f(Js) 
f'(Js) 



log uj + 2\J 'A logo; 
2 log 3/2 u 
A 1 / 2 ' 



(5.18) 



where indeed we can see that at J s we have /" > 0. Now the invariant amplitude 
becomes 

>■ dX 2 

n=0 J Vo 



i + 



7rA(X 2 ) 1 / 2 e log u; g 2 y/A(\*)logLj 



r 3 / 2 , 



27rVlog 3/2 o; M ^ 



4(A 2 ) \| log"" - ,C 
logo; 

Q 2 rfA 2 



-a s (A 2 )A(A 2 ) 1 / 4 - 



1 + 



logo; 



oo 

•e 2 ^ (A2)log - E" B ^&,n(«.(5),l)(?/||O a ' )(J ' B) lb>( w ,) • (5-19) 



n=0 
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Because the saddle point fixes the value J — n = J s — n — 1+ y j^j, which fixes the 
number of internal derivatives inside the gluon operator, recall (|3.18|) , it is clear that 
the reduced matrix element in ( |5.19| ) is the conventional (forward) gluon distribution 
function in moment space with a shifted moment label J s — n. Despite the formal 
summation over n, which is essentially the only difference introduced in this limit by 
the non-forwardness, only one non-perturbative input is needed, which is the gluon 
density. This means that although in general there will be complicated dependence 
on v in the double distributions, (in high energy and hard scattering limit) it is 
nontheless perturbative. 

At leading order the n summation collapses when we take the lowest order Wilson 
coefficients given in ( |4.60| ) with the flavor averaged quark charge rif < e 2 >= J2f e q, 
which actually could be done even before taking the saddle point approximation. We 
went to the saddle point first in order to make manifest the fact that there would be 
only one soft input, the forward gluon density, even with the formal n summation. 
Thus, subsitute the result from equation |4.60| we have 

. f Q2 ^ as{X 2 )A(X 2 ) l/4 ■ 1 e 2 ^^. (5.20) 

Jul X 1 _j_ / A(X 2 ) 

1 + V 

Because of the pinching of J s — n, the final evaluation of ( |5.20| ) is the same in both 
forward and non-forward cases, in the high energy limit under leading logarithmic 
approximation. Define the integral 

<2 2 d\ 2 ,,,, 1 



a s (X 2 )A(X 2 )^ :== eV^ 2 ) 1 ^- , (5.21) 



M5 A 2 1 

1 + V log^ 



with A{X ) given in equation 5.10 , which is not very big even for very hard processes 
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due to the double logarithm. Therefore in the high energy limit where u is very big, 
we have 

1 . (5.22) 



1 + 

We make the change of variable 



\oguj 



log(A 2 /A 2 ) = log(Q 2 /A 2 ) - u , (5.23) 

which implies that 

u = log(g 2 /A 2 ) - log(A 2 /A 2 ) = log(g 2 /A 2 ) . (5.24) 

It is clear that u > and is not very large either. To simplify the notation, let 

AC A 



a = 



bo 

b = log(g 2 /A 2 ) 

c = l0 >g(Q7A 2 A 
g (,log(^/A 2 ) ) ' 

where we make the note that a and b can be seen as moderate compared with the 
double logarithms u and c, which are not big. Then 

— 2, . ox 

U 



log(^/A 2 )/ s ^logM/A 2 )^ log(Q 2 /A 2 ) 
= c + log(l--) 

And the integral / becomes 

2 

Jo bo(b-u) Abe 

= du A acf/Hl - — ) e 2 v / ^^a-^) 

Jo bo{b-u) Abe 



fac) 1 / 4 
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where the second integral II is defined as 

2 

// = d„-J-(l - -£-) e-f . (5.27) 

Jo o — u 4oc 

As we have discussed earlier, because the relative size of a and c and that logcu is 



considered big, y log ^ a is actually a big number. 
Rewrite (f[27|) as 



1 f 10 ^ , 1,1,. u, . I 

J/ 



-(— (1 - -) + 1 e^v —r- 

l-£ v 4c b J Ac 

b 



1 flog^ 2 



and let 



we have 



« = iog(g7A 2 

b log(g 2 /A2) 



g = ; = ^Z/,: j , (5.29) 



II = / o d 2 ( f -VSS + (i - J • (5-30) 

The z integral can be evaluated using the formula of special integral exponential 
function IEi(z) 

/ e- Xz = e-\lE t (\)-IE l (\-a\)), (5.31) 

JO 1 — 2 

where the IEi(z) function has a branch cut discontinuity at the negative real axis in 
the z-plane (i.e., cut at z G (— oo,0]), and has the following properties: 

IEi(0) = -oo, 
7^(0.3725) = 0, 



Ill 



and while 



d z IEi(z) = -, (5.32) 

z 



we have 

IEAz) 

— — « 1 5.33 

e z 

when z > 0(1). 



Therefore the second term in equation |5.30| does not contribute in the leading loga- 
rithmic approximation, and we have 



1 I / --A^^I^M' 

II = — , 1 - e V 



4c ylogua/c 
1 

4y/\oguac ' 



where in the last step we have used the fact that J log uia/c is a big number. 



(5.34) 



Subsititute these results into the expression of the invariant amplitude ( p.2U| ) we 
have 



/ * (ac)V 4 , 1 



V log 3/2 uj b Ay/logcoac 
Since ac = ^A(Q ), the final expression of the double distribution T 2 is 

P-Q pr~?T (7? ^ in f <e l> 
— Jl-(T 2 {Q 



. (5.35) 



.y||dW^||p>^a, C V^^. (5.36) 
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In the very high energy and hard scattering limit, the asymmetric reduced matrix 
element becomes the forward gluon density 



( P '\\d^ J ^\\p) -> sji^c (p\\d i9)iJ ^\\p) 

ex /K(x#)), (5.37) 

with x = —. We can see clearly that the leading high energy (small-x) behavior is 
exactly the same as that is given by a forward direction leading logarithmic analysis 



(see section p..5.3|) . 
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Chapter 6 



Conclusions and Open Problems 



We have formulated a general operator product expansion of a non-forward and un- 
equal mass virtual Compton scattering amplitude. We have found that, because of 
the non zero momentum transfer, the expansion now should be done in double mo- 



ments with respect to the moment variables defined in equation |3.13| . The double 
moments can be parametrized as products of the Wilson coefficients, which can be 
computed perturbatively at the hard scattering scale, and the non-forward matrix 
elements of new sets of quark and gluon operators, namely, the double distribution 
functions. These operators in general have total derivatives and they mix among 
themselves under renormalization. They obey a well-defined renormalization group 
equation that can be solved formally. We have calculated the (equivalent to) evolu- 
tion kernels of the double distribution from which the anomalous dimension matrix 
of these operators can be extracted. 

In the high energy limit, we have found the leading contributions of the anomalous 
dimension and used them to solve the resulting renormalization group equation. We 
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have recovered the conventional double leading logarithmic (DLL) analysis results and 
have made the connection that in fact the double distributions are proportional to 
the conventional forward gluon density in this limit. This justifies previous analyses 
using DLL approximation and forward gluon density on non-forward processes like 
the exclusive diffractive vector meson production (e.g. |I7|). Furthermore, we have 
developed the formalism to proceed in principle beyond DLL to next to leading order 
and/or relaxing the leading logarithmic approximations. 

Our analysis is quite general. We can go from deeply inelastic scattering (DIS) 
(the forward case, corresponding to v — 0) to deeply virtual Compton scattering 
(DVCS) (corresponding to v — 1) in principle by continuously varying the moment 
variable v. However, the analyticity/analytical continuation in z/, as well as the 
relevant dispersion relationships, still need further investigation. At the same time, 
diffractive vector meson production {y close to 1) can also be studied by replacing 
the second quark-photon vertex with effectively a (light cone) wave function of the 
vector meson |17|, [31]j . Also, here because the time-like four-momentum of the vector 
meson, further study is needed in extending the discussion to the time-like region. 

There are also other details of the analysis still need to be filled in. The ex- 
plicit form of the anomalous dimension matrix is not yet written down, although 
as mentioned earlier, a change of base to a linear combination of the quark and 
gluon operators |3.11| using Gegenbauer polynomial as coefficient might simplify the 
situation. The dispersion relationship between the invariant amplitudes Tj and the 
double structure functions need also to be clarified. For example, the exact physi- 
cal meaning of the fact that different k_ poles are being picked out at different k + 
integration regions. The relationship between this and the transition/interpolation 
between Altarelli-Parisi evolution and Brodsky-Lepage evolution is also an important 
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topic to study in depth. 

More experimental data are to come from DESY, Jefferson Lab, SLAC, CERN 
and RHIC, the study of the physics of non-forward processes will continue to be an 
interesting and important field. 
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Appendix A 

Quark-Gluon Anomalous 
Dimension 



In this appendix we show the details of the calculation that leads to our form of the 



quark-gluon anomalous dimension in equation [4.53 



<qg 



In the LC gauge, the value of the quark-gluon diagram, shown in figure |20 

d 4 k 

K^ab^ba) / 7 

a, b 



IS 



It- / 



$ $1 $ J 



(a.i; 



where O q is the quark vertex from ( |4.23| ) and there is an extra factor of (—1) coming 
from the fermion loop. After the color algebra we obtain 

,1. f d 4 k (2k-r) J ^ 1 rlB a p( y k,r) 



(2vr) 4 (fc 2 + ie)((Jfei-Jfe) 2 + ie)((k-r) 2 + ie) ' 



(A.2) 



where we have defined 



B af3 (k,r) = Tr( la (ft- /) 7+ ^(fa- ft)) 



(A.3) 
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Almost identical to the calculation of the gluon-quark section, after the fc_ integral 
we have 

ytf _ <x*l fdk 2 ( r + dk + k + M (1) 

4) 



where .B^ and £>( 2 ) are B a/3 (k,r) evaluated at the poles fcL 1 ^ and /ci 2 " 1 , respectively, 
and again with 

k 2 k 2 

k (i) = J_ = & A; 2 = 2A; + A;_ - A; 2 = 0, 

k r2: . ~ . — ^ & A; 2 = 



2(fc+-fci+) 2Jfei+(l-w) 
We have 

= Tr( 7a JWi" 0) - Tr( 7a / 7+ JM ft- 0) 

= #J,a/3 — B II:Ct p . (A. 5) 

Same as in the main text, we will pick up only the terms that are logarithmic 
divergent in the transverse momentum (A; 2 terms) and thus we need only look at the 
transverse components (a/3 — i,j — 1,2 terms). 

For the first term, we have 

Bl,af3 = TV (7 a /^7+ #7/j(#L- #)) 

= Tr( 7a (7 + ,^-7 + JMfr" fl) 

= 2A;+Tr( 7a # 7/J (fr- #)) - k 2 Tr{ lal+lp {fa- #)) 

= Sk+ik^h - k)p + k p {k 1 - k) a - g af3 k ■ (k x - k)) 

-A; 2 Tr( 7 _7 + 7_7 + (A; 1 - fc)_) - k 2 Tr( 11+11 4h - fc)+) , (A.6) 
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where we have used the identity of gamma matrices M 

Tr{4 PVP)= A{{A ■ B){C ■ D) — [B ■ D){A -C) + (A-D)(B- C)) . (A.7) 

At the pole fc_ where k 2 = 0, we have, for the ij component, 

B^ = 8k + (-kikj - kjki - Qijk ■ h) 

k 2 

= 8k + (gijk 2 - g t ~) 

= 8k +9i:j k 2 (i - i-) 

= Ak 1+9ij k 2 {2uj - 1) , (A.8) 

while at the pole k_ , where factors of transverse momentum k can come from either 
k 2 or kikj terms, we have 

B i% = 8k + (-k i k j -k j k i -g ij k-(k 1 -k)) + Ak 2 g ij (k 1 -k) + 

1 k 2 Ak 2 
= 8/c+(^fc 2 - -g ijT =— ) - - w) 

Z 1 — U) 1 — LU 

= Ak 1+gij k 2 (2uj - — 1) 

1 — LO 

= -Ak 1+ g l3 k 2 U}2 \ {l ~ . (A.9) 

1 — UJ 

On the other hand, for the second term we have, by using the anti-commutation 
relationships of the Dirac matrices ( |1 . 1 1 7| ) , 

Bn, a p = Tr(7 Q jhf+ ^75(^1- #)) 

= 2 7a Tr( 7+ # T/ j(ft- ¥)) -Tr(/ TaT+ #7/j(ft- #)) 
= 2 7a Tr( 7+ ^(ft- 0) " 2<? a+ Tr(/ # 7/9 (fr- #)) 
+Tr(/ 7+7a ^(fr- m 
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2 7a Tr( 7+ m - 2g a+ Tr(f # T/9 (fr- )6)) 

+2k a Tr{h+lp{h- }6))-2g a/3 Tr(f 1+ #)) 

+2(h - k) a Tr(f 1+ fag) - Tr(/ 7+ # 7/9 (ft- #) 7a ) , (A.10) 



which gives 



■8/1,00 = J a Tr(j + JijpiJti- #)) ~ 9a+Tr(jf # 7/3 (#i- #)) 
+A; a Tr(/ 7+7/3 (# 1 - #)) - <7 a/3 Tr(/ 7+ #)) 
+ (fci-fc) a Tr(/ 7+ # 7/3 ). (A.li; 



At the pole fci 1 ^, the first term of ( A.11J) , after taking only the logarithmic divergent 



contribution and only the transverse component, becomes 

7a Tr( 7+ ( 7 _A; + - 7 • fc^fr+fa-A;)- + j-(h-k) + - 7 • (fcj-*))) ^ ; (A.12) 

the second term is always zero because gj + = 0; the third term now becomes 

4£; a (r + (£;i-£;) / 3 + gp + r ■ (fa-k) - rp(ki-k)+ -Akik 3 r + = 2g ij k 2 r + ; (A.13) 

the fourth term is now 

—4gij(r + k ■ {k\ — k) + r ■ {ki — k)k + — k ■ r(ki—k) + ) 

k 2 k 2 k 2 
= -^(r + ^-r + ^-^-r + (h-k) + ) 

= -2g ij r + k 2 (±-l-—) = 0; (A14) 

UJ UJ 

and the last term is 

- ^Tr(/ 7+ f l3 ) = -Ahir+kj) = 2 9ij k 2 r + . (A.15) 
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Thus we obtain 

B% = ^k 2 r + = Ak 1+9lJ k 2 u . (A.16) 

(2) 

At the pole k_ , the first and second terms are both zero for the same reason; the 
third and fifth terms will take the same value as they do at the fc_ pole; for the 
fourth term we now have 

—4gij(r + k ■ (ki — k) + r ■ {k\ — k)k + — k ■ r(ki — k) + ) 

h 2 h 2 r h h 2 h 2 

All — — \ ' + + / UL \ I I t \ 

= - 4 ^ (r+( "2(T^) + (T=^ + " r+( ~^-(I^) )( 1 )+ 

= -2^(1^+1^+1) 

= -A 9ij k 2 r + -^— . (A.17) 

1 — LU 

And we obtain 

>(2) 



B n\ij = ^k 1+ g i:j k v - Aki+gijk - 



to 

2 



= -ik 1+gij k z - . (A. 18) 

Collecting terms we have 

fig' = B£?.-B$y=4* 1+9e i 2 (2w-l-i/) 

B%> = - 3%, = -ik 1+9ij e> ^-^-"v . (A.19) 



Substitute this back into ( A.4j) we have 



*i+ ^ + (/c 1 -/c) +(2fc _^^„ v „^ 2 + (l-^) 2 -u;z/ N 



<ifc 2 / f r + dk+kj 



(ki-r) + v /+ + l-u 

— ~2g ii kLv n / -=-( / 2w - z/ V 7 "^ 1 (2w-l- 

2tt2 yy 1+ J k 2 K Jo v K ' K 



1 ^_{2uj-u) m (uj 2 + (1-uj) 2 -uju) \ . (A.20) 
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After factoring out the logarithmic divergence in fc 2 we again recover the result from 



(EH) 



'19 



2tt2 



2g ij k J 1+ v n J^du{2uu-v) M (-Q{v-u)^uj(2u- 1 -u) 

+®(u]-u) — (lu 2 + (1-luY-luu)\ . (A.2i; 
1 — f / 
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Appendix B 

Gluon-Gluon Anomalous 
Dimension 



The computation of the gluonic sector follows similarly. In this appendix we will show 
in some detail the calculation of the gluon-gluon anomalous dimension. 

Using standard Feynman rules and light-cone (LC) gauge we can write down the 
value of the gluon-gluon triangle diagram (figure ETf) as 



(2tt) 4 (k 2 + ie^h-kf + te)((k-r) 2 + ie) 
■D^{k - r) Offish, k-r) D a , a (k) (B.V 



where the gluon vertex 0^ J,n - ) is again taken from ( [4.23 ) and T denotes the triple-gluon 
vertex. 

Performing the color algebra with 

^Z(fijkfik'j) = —CA5k,k' (B-2) 

ij 
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and setting the index k = k' we arrive at 

d^k 2(2k - r) J ^- 2 r\Al v (k, r) 



7 



y A J (2ir) 4 (k 2 + ie)((h-k) 2 + ie)((k-r) 2 + ie) 1 ' ' 



where 

A% = T, aS D, s (h - k) TsvpDwik - r) Vjj, a ,(k, k-r) D a , a {k) (B.4) 

with 

and similar Dppt{h — r) and D a i a (k). Again the LC null four vector n is defined with 
properties n 2 = and n ■ v = v + for any four vector v. 

The one point we need to pay extra attention to is the form of the gluon vertex in 
( |4.23| ) and thus in (|B.4|) . The full tensorial structure of the LC gluon vertex in the 



non- forward case, generalized from that in pi, is 



Vp, a ,(k,k—r) = g^am-k^k—^—n^k—rJk'pn^—n-kn'^k—ry^k^k—^n^n^, (B.6) 

and it is to be contracted with gluon lines D a r a (k) and D^ik — r). 
We are now going to show explicitly that 

Dppik - r) V$, a ,{k, k-r) D a , a (k) = V° a (k, k-r) 
Dpf3<{k - r) gp a ,n-kn-{k - r) D a , a (k) = V$ a (k, k - r), (B.7) 



thus we can substitute the full gluon vertex by its first term because of the LC 
projectors of the two gluon lines connected to the vertex and use ( 4.23 ) rather than 
the full expression in (|B.4|) . 
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We split Vp, a ,(k, k — r) into two parts, with the gpi a i term being one part and the 
remaining three terms being the other. That is, we write 

Df,p(k - r) V° a ,(k, k - r) D a , a {k) = V™ - vfc 2) , (B.8) 

where 

V]f a l) = Dpp(k-r)gp> a >n-kn-(k-r)D a i a (k), 

v pa 2) = D P/3'(k - r ) ( n '(k - r )kf3>n a > + n-knp>(k - r) a > -k-(k- r)np>n a i) D a , a (k) . 
We have 

= {n-kn-(k-r)gp a , - n-k(np(k-r) al + (k-r)fin^)) fa a - — ka ^ ka Wq ^ 
= n-kn-(k-r)gp a - n-k{np(k-r) a + {k-r)pn a ) - n-(k-r)(npk a + kpria) 

+ n-(k — r)npk a + k-(k — r)nfsn a + n 2 (k—r)fsk a + n-k(k — r)pn a 
= gp a n-kn-{k — r) — n-(k — r)kpn a — n-knp{k — r) a + k-(k — r)npn a 
= V p 9 a (k,k-r), (B.9) 



and 



= fa, - Mk ~ r ^-^) r) ^ ) {n - {k - r)k ^' + n - kn ^ k - r) «' 

-k-(k-r)np>n a >) D a > a (k) 
= {n-{k—r)kpn a i + n-knp(k — r) a i — 2k-(k — r)n/3n a ' — n-k{k—r)pn a i 
'n-(k-r)np(k-r) a , + n 2 (k-r)p(k-r) al ) 



n-(k — r) 
k-(k-r) 
n-(k-r) 



[n-{k-r)npn a i + n 2 {k-r)pn a ,)) D a , a (k) 
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(n-(k — r)kpn a > + n-knp(k—r) a > — k-(k—r)npn a i — k-(k—r)npn a i 
—n-k(k—r)pn a i — n-knp(k—r) a > + k-(k — r)npn a r) D a > Q (k) 

(n-(k—r)kpn a i — 2k-(k — r)npn a > — n-k{k—r)pn a i + k-{k—r)npn a i) 

(fi a 'k a ~\~ k a in a \ 
9a ' a —k ) 

(n-(k—r)kpn a — 2k-(k — r)npn a — n-k(k — r)pn a + k-(k — r)npn a 



n kpk a + n-kkpka) H (n /cgfco, + n-kkpk a ) 



n-k n-k 
+ (n 2 (k — r)pk a + n-k(k — r)pn a ) ^— - — -(n 2 npn a + n-knpn a )) 

Tl - K 

= n-(k—r)kpn a — k-(k—r)npn a — n-k{k — r)pn a — n-(k—r)kpn a 

+2k ■ (k —r)npn a + n-k(k — r)pn a — k-(k—r)npn a 
= . (B.10) 

It is now clear that we have established (P-7| ), and explicitly we have 

Dpp>(k - r)V$, al (k, k-r) D a , a (k) = V^Jk, k - r) 

= gp a n-kn-(k — r) — n-(k — r)kpn a — n-knpik — r) a + k- (k — r)npn a . 

(B.H) 



The other factors in the definition of A 9 (equation B.4|) are the triple-gluon 
vertices involved in the diagram: 

F-ya^ = g-y a (2k — ki)^ — g a ^{k\ + k)^ + g^(2ki — k) a = — T^^ a , 

Fsvp = gs v {^ki-k—r)p — g v p{2r—k 1 -k)s + gps{2k-kx—r) v = —Yp v s. (B.12) 



Thus we can rewrite ( B.4j) as, since = -D57, 

A % = ^aftD^ s {ki-k)T S ui3Dpi3'(k-r)V^ a ,(k i k-r)D a f a (k) 
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= T ia5 D l5 (ki-k)T 5l/l3 V$ a (k,k-r) 

= {—^5vp)Ds 1 {—Y ia ^)Vp a = Yp v sDs 1 V llia yp a . (B.13) 
Expanding the above we obtain the full expression of the "gluon triangle" as 

A 9 ^ = r 'f3v5Ds 7 r JIMX Vp a 

= n-kn-(k-r)[((k-2k 1 )-(k-2k 1 +r) + (fci+Jfe) -(fci+fc— 2r) ) 
+ (fci+A;-2r) /1 (A;-2A;i) I/ + (k-2k 1 +r)^(k 1 + k) u 
+ (fc-2fci) M (A;i-2A;+r) 1/ + (fa+fy^fa^k+r^ 

+ (k 1 -2k)^k 1 + k-2r) u + (fci-2fc) M (A;-2A;i+r) 1/ + 4(^-2^)^(^-2^+^^] 

- n ' kn ( '} k ~^) [((fci + fc)-(fci-fc)w(fci + fc-2r) + (fci-fc)-(fci + fc-2r)7i-(fci + fc))^ 
+ n-(A;i+A;-2r)((A;i-A;) A1 (A;-2A;i) I/ + (A;i-2A;yA;i-A;) I/ ) 

+ n-(fci + fc) ((/ci-A;) M (A;i-2A;+r)^ + (fe — 2/c x +r) M (fci — /c)„ ) 

+ n-(k-2k l )(k 1 -k) fl (k l -2k+r) u + n-(k-2k 1 +r)(ki-2k) fl (k 1 -k) u 

+ 2n-(k 1 -k)(k 1 -2k)^(k 1 -2k+r) l/ 

+ ((k 1 -k)-(k 1 +k) + (k 1 -k)-(k-2k 1 ))n lt {k 1 -2k+r) v 

+ (k 1 -k)-(k 1 + k-2r)n fM (k-2k 1 ) v + (k-2k 1 )-(k-2k 1 +r)n fl (k 1 -k) v 

+ ( + A; -2r) + (A?i — fc) • (fc — 2/ci+r) ) (k^k)^ 

+ (k-2k 1 )-(k-2k 1 +r)(k 1 -k)^n I/ + (k 1 -k)-(k 1 + k)(k-2k 1 +r%n l/ ] 

— n-(k-r) [n-(k — 2ki) k-{k — 2ki+r)g^ u 

+ n-{k-2k 1 ){k„(k 1 -2k+r) u + (k 1 + k-2r%k u ) 

+ n-(k 1 +k-2r)(k 1 -2k) li k v + n-k(ki-2k) IM (k 1 -2k+r) v 

+ (k 1 +k)-(k 1 +k-2r)n lx k l/ +k-(k-2k 1 +r)n lx (k 1 +k) l/ +k-(k 1 +k)n lx (k 1 -2k+r) u 
+ k-{k-2k l +r){k 1 -2k) fx n u ] 
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+ —71 7T [n-(k 1 + k-2r) ( n-(fc-2fc 1 )(fc 1 -fc) /i fc I/ + n-(k 1 -k)(k 1 -2k) fl k l/ ) 

+ ( n-(fc-2fci)(fci-fc)-(fci + fc-2r) + n-(fci + fc-2r)(fci-fc)-(fci + fc) 

+ n-(fci+A;)(A;i-A;)-(A;i+A;-2r) ) 
+ k-(k-2k 1 +r) ( (ki-k)-(ki + k) n^n u + + +n-(A;-2A; 1 ) )n^(ki-k) u 

+ n-(k — 2ki)(ki — k)^n u + n-(ki — k)(ki — 2k)^n u ) 
+ ( (n-(k-2k 1 )k-(k 1 -k) + n-k(ki-k)-(k 1 +k) + n-(k 1 +k)k-(ki-k))n ll 
+ n-k(n-(k-2k 1 )(k 1 -k) li + n-(k 1 -k)(k 1 -2k) li ) ) (fci-2fc+r) v ] 
-n-A; [ n-(k-2kx+r) (fc-2fci)-(A;-r) 57^ + (k~2k 1 )-(k-r)n II (k l -2k+r) l/ 
+ n-(k-2k 1 +r) ( (A; 1 -2A;) M (A;-r) l/ + (k-rj^h+k),,) 
+ {n-{k l + k) (k-r% + n-(k-r) (h-2k)^ ) (fci — 2fc+r) v 
+ ( (A;-2A;i)-(A;-r) (fci+fc-27% + (fci+A;)-(A;i+A;-2r) (fc-r)„ 
+ (fc-r)-(fci+A;-2r) (fci-2fc) M ) ] 
+ 7J~~ i \ [ ( (fc-2fci)-(A;-r)(n-(A;i+A;-2r) (A*-*;),* + (A; 1 -A;)-(A; 1 +A;-2r) n M ) 
+ (n-(fci+A;-2r) (fci — A;) • (fci+fc) + ra-(fci+A;) (ki-k)-(k 1 + k-2r) ) 
+ (n-(fci+fc-2r) (fci-A;)-(A;-r)+n-(A;-r) (fci-fc)-(A;i+A;-2r) ) (fci-2fc) M ) ra„ 
+ n-(k-2k 1 +r) ( (n-(fci+A;) (fc-r) M + n-(fc-r) (fci-2fc)^ 

+ (fc-2fci)-(A;-r)n M )(A;i-A;) I/ 
+ ( (A;-2A;i).(A;-r) (h-k)^ + (h-^^h+k) (fc-r)„ 

+ (k 1 -k)-(k-r) (fci-2fc) M )n I/ ) 
+ n-(fci-A;)(n-(A;i+A;) (A;-r) M + n-(k-r) (fci-2*;)^ 
+ (fc-2fci)-(A;-r) n M ) (fci-2fc+r)„ ] 
+ k-(k— r) [ n-(k — 2ki) n-(k — 2ki+r) g^ v + (fci+A;) • (fci+fc — 2r) n M n^ 
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+ n-(k-2k 1 + r) n^kx+k),, + (n-(k-2k 1 ) + n-(h+k))n lx (ki-2k+r) I/ 

+ n-(k-2k 1 ) (k 1 + k-2r) f ,n„ + (n-(k 1 + k-2r) + n-(k~2k 1 +r) ) (ki-2k)^n v } 

— ~~~7~r~~J\ [ {n-{k-2kx) +n-(k l +k))n-(k-2ki+r)nJk 1 -k) u 
n-{k\—k) 

+ (n-(k-2ki) + n-(ki+k) ) n-{k\ — k) n fJ ,(k 1 -2k+r) u 

+ (n- (k!+k-2r)+n- (k-2h+r)) (n-(k-2ki) (ki — k)n+n-(ki—k) (ki-2k)^)n v 
+ ( (n-(fc-2fci) +n-(k 1 + k) ) {k 1 -k)-{k 1 + k-2r) 
+ (n-{ki + k-2r) + n-(k-2k 1 +r) ) (k 1 -k)-(k 1 + k) ) n^n v } , (B.14) 

where we have used n-n = n 2 = 0. 

Similar to the calculation of 7 gg we perform the k- integral in ( |B.3|) first. The 
discussion of different regions of k + and the corresponding pole values goes almost 
identical. We have 



. 2 ,, i d±k 2(2k- r y^rlA^(k,r) 



(2tt) 4 (k 2 + ie)({ki-k) 2 + ie)((k-r) 2 + ie) 



(2tt) 2 J r \ r+ k 1+ Jo 



+ , ^ ; r + dk + (kt - k) + A^(2k - r)J"V 



- r) + 



OsCa [d£( f r+ dk+k+ M (1) 

2tt 7 fc 4 Wo r + A; 1+ V ; + + ^ 



+ 



fci 



+ dk + (k\ — k) + 
r+ -h+ih - r)+ 



(2k - r) J ^\lA^A 



7^+7^, (B.15) 



where A 9 '^ and A 9 '( 2 ) are A^(fc,r) evaluated at the poles fc_ and k^\ respectively, 
and again we have 
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k 



(i) 



2k, 



k 



(2) 



k 



-k' 



& k 2 = 2k + k. - k 2 = 0, 



& k 2 



;b.i6) 



2(Jfe + -Jfei+) 2/ci + (l— 

To extract the gluon-gluon anomalous dimension, we do not need to calculate 
further the full expression B.14 . Rather, the projection of A 9 onto the transverse 
directions will suffice. That is, by effectively treating the gluon vertex to be 



Vfj = QijU -kn-{k — r) 



;b.i7) 



we only need to keep in equation |B.14| those terms proportional to gij and kikj 
— \9ijh 2 - Thus we define 



A\ '■ = g^Ai + kikjAn 



(BA8) 



and calculate Ai and An separately. Note as stated in the main text, we may choose 
k\ and r such that only their plus components are non-vanishing, that is, 



fci = (*;+, 0,0), r = (r + , 0,0) 



(B.19) 



which also leads to 



k\ ■ r = k\ 



Reading off (and simplifying) from the full expression of A 9 we arrive at 



(B.20) 



Aj = n-kn-(k-r)(2k 2 -2k 1 -k-k-r) 

+ v ,/ (2k 2 n-(k 1 + k) - 2k 2 n-r - 2n-(k l + k)k-r) 

n-(ki-k) 

-n-(k-r)n-(k-2k 1 )(k 2 -2k 1 -k + k-r) -n-kn-(k-2k 1 )(k 2 -2k 1 -k) 



— n- k{n-r{k 2 — 2k\- k — k ■ r) — n-{k — 2k\)k-r) 
+ ( y k 2 -k-r)((n-( y k-2k 1 )) 2 + n-r{k-2k x )) . 



(B.21) 
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Same as the discussion of the quark-quark anomalous dimension, we are only inter- 
ested in the leading powers of k, which is also quadratic, and, from the pole values, 
can only come from terms proportional to either k 2 or k_. Keeping only such terms 
we obtain 

Aj = k 2 Aj + {k 2 - 2ki ■ k)A] + r+k^Aj , (B.22) 

where 

A a j = 2\k\ — 2k 1+ k + + 2k'f + - k l+ r+ + 



(ki - k). 

, f UJ 2 + UJ 3 — ujv — 2uj 2 v + LOU 2 \ 

= 2k 2 1+ iuj 2 -2uj + 2-u + — 1 

A b j = Ak 1+ k + -k\- 2k 1+ r + - k + r + = k 2 1+ (Auj - uj 2 - 2v 
Aj = -2 (jfcj. - k + r + - 2k l+ k + + 2k 2 1+ + ~ k >+ - k ^ k ^ 



UJV) 



(h - k)+ J 

,9/9 UJ 3 + UJ 2 — UJ 2 U — UJV s . 

-2k 2 (uj 2 - ujv - 2uj + 2 + ) . B.23) 

1 — UJ 



At the k_ pole, from (|B.16|) we have 

,2 



k 2 = 0, 2k v k = —, (B.24) 

UJ 



thus 



uj 2uj 



k 2 

-— k\ , (Auj — uj 2 — 2v — ujv) 

UJ 

vk 2 , 9 , 9 u; 3 + u; 2 — cuV — ujv . 
-— A; 2 , (cj 2 - wi/ - 2cj + 2 H ) 

UJ 1 — UJ 

-— /cf , (4u; - cj 2 + viuj 2 - 3uj - ujv H -)) .(B.25) 

uj 1 — UJ 



At the fci 2 '* pole, on the other hand, from ( P-16Q we have 



k 2 

k 2 = 2k\-k = — — = — , (B.26) 

1 — to 1 
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and thus 



A {2) _ k 2 uk 2 

= -- k( + (2(u-2u + 2-u^ ) 

1 — UJ 1 — UJ 

. 9 UJ 3 + uj 2 — u 2 u — UJU. . 

-u(io 2 - iou - 2oo + 2 H )) 

1 — a? 

^ 2 „,2 / 2 „ „ „ tu 3 + ^ 2 - 2uis(2u - I/) . „„s 

= — — 2&; 2 , (cj 2 -2a; + 2 - 2i/ H -) . (B.27) 

1-a; 1+v l-u 



On the other hand, from the full expression of A 9 ( |B.14| ) we have 



An — 10n-kn-(k—r) — 2n-kn-(k—r) — n-(k—r)n-(k+4r) — n-kn-(k—5r) 

n-kn-(k-r), , n . ,, , ,, 

+ /; (w-(fci-fc-2r) + ra-(fci-fc + 2r)) 

n-{k\—k) 

= 8k 2 + - 8k + r + + 4r 2 + 

= Ak\ + (2uj 2 - 2ujv + v 2 ) . (B.28) 

Note that A// has the same value at both poles. 
So now we have 

■A-ij = QijAi + kikjAjj 
= g^Aj -h 2 A H ) 

= gij(Aj - k 2 k 2 1+ 2(2uj 2 - 2uu + z/ 2 )) . (B.29) 

Substitute the above result with the values of Aj at different poles into ( B.15|) , we 
have 



lj(1) a s C A fdk 2 r+ dk + k+ M nA g,(i) 

a s C A fdk 2 r^,M,„ vr^n _ , „m 1 



2tt 
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and 



2tt J k m 1+ Jo V ; 2u 



, 1. 9 . 9 cu 3 + u; 2 — viuj 2 + cj) . . 
• Uuj - uj 2 + v[uj 2 -3U-UU+ i '-)) 

UJ 1 — UJ 

-2(2u 2 - 2cjz/ + v 2 )) 
oi s Ca f dk 2 



1 (l, , o , u; 3 + cj 2 — 2u; 2 z/ . 9 N \ 

- -(4u; 3 -u; 2 + 4u;) + — 3 cj 2 + B.30 

2 \ u I — uj 



1,(2) _OsCa f d£_ dk + {k 1 -k) + _ M g , (2 ) 

- 27T J k A Jr + k 1+ (h-r)+ [ )+ + ij 



^ / f / d -^k{- 2 (2. - uY^ 9iA Af - h 2 A n ) 



asCA f^xJ'd^-uV^ l ~ u 

J k Jv 



2tt J k 1 ^"^Ju v ' 2(l-i/) 

2 , 2 o n o tu 3 + tu 2 - 2u;z/(2cj - i/) 
(cj 2 - 2u; + 2 - 2i/ + 



1 — UJ 1 — UJ 

2 o ,.2\ 



-2(2uj 2 - 2uju + u 2 )) 

cx s Ca f dk 2 
2vr 



\ J- j n— 2 ^jn 



f dk r 1 
J -j^ r 2g ij k( + J du{2u-v) 

1 ~/2(l - + cu 2 ) 2 + + uj 2 ) - 2(1 + uj 3 )Y 



\ — u\ I — CO 



(B.31) 



Same as in the case of the quark-quark anomalous dimension, the uj integration 
is divergent at the end point uj = 1. This divergence is again cancelled by self-energy 
graphs. However, in this case for each gluon propagator connecting to the vertex, 
there are two of them, a gluon loop and a quark loop (figure |22|). Their values in the 
LC gauge can again be taken readily from j|7| as 

Z g {x) = l + ^2C A -(-/o-log|x' " 



2tt e 01 1 12' 

^2C A \- n -l- 
2tt e U 



Z q {x) = l + £.2C A -(--£). (B.32) 
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Again the self-energies depend on the longitudinal momentum fraction x, with Iq the 



LC colinear divergence (see equation f4.51|). Using the same argument as that was 



used for the quark-quark graph, we need add to the gluon-gluon anomalous dimension 
the following contribution from self-energy diagrams: 

\{z g {\) + z q (i) + z B (i-u) + z q (i - u)) = ^C A ~(-i - i i g(i - 1/) + 11 - . 

(B.33) 

Identifying - as the logarithmic divergence in transverse momentum and that 

11 - 5/ = J-(ii - ^ ) s _^ 6o (B . 34) 

12 18 4iV c v 3 /; 4iV c v ; 

with 6o the leading coefficient of the QCD f3— function, we can see that we have shown 
that the gluon-gluon anomalous dimension is indeed given by the same expression of 
j gg as in Q4.53| ), that is, 



7* = ^20«A£i/ n j(W^ 

+ i 3(a; + — + i/ 

1 — UJ \—V 1—UJ 1 — UJ 

+ (l-u)5(l-uj)(-2I -\og(l-u) + ^)\ . (B.35) 
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Appendix C 

Gluon-Quark Anomalous 
Dimension 



In this appendix we show the details of computing the gluon-quark anomalous dimen- 
sion. The gluon-quark transition comes from the contribution of the second triangle 
diagram, as shown in detail in figure |19[ By the same Feynman rules and in the LC 
gauge, we can write the value of the gluon-quark diagram as 

. n (J,n)(+) (h h , -iP a »{k) 

^ k r > k 2 + i£ 

d 4 k 2(2k - rY^rlB^k, r) 



19 Cp5aa ' ' (2vr) 4 (k 2 + te^ih-k) 2 + ie){{k-r) 2 + ie) ' ^ 



Similar to the calculation of the quark-quark diagram (see section 4.2. 1| ) we have set 



the color label of the two quark lines to be equal (a = a'), and defined 

B"(k,r) = 7„(ft- tehuD u p(k-r)V 9 a D a ^k). (C.2) 
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The LC gluon projectors are 

D v p(k - r) = g vP 
D a ^(k) = g a ^ 



n u {k - r) p + (k - r) v np 
n ■ (k — r) 

n ■ k 



(C.3) 



From our convention of the gluon vertex ( 4.23|) , by exactly the same argument as 



that for the gluon-gluon diagram (see the derivation of equation |B.7| in appendix B) 
we have 

D uP {k - r) V$ a (k, k - r) D a ,(k) = V^(k, k - r) 

= g Vfl n-kn-(k — r) — n-(k — r)k v n^ — n-n u (k — r) M + k- (k — r)n u n^ . 

(C.4) 



By the same steps we followed before in computing the k_ integral and picking 
up different poles at different k + regions, we have 



7, 



99 



a s C F f dk 2 ( r+ dk + k 



2tt 



k \Jo r + k 1+ 

dk+(h - k) 



+ 



(2k - r y^ r n Bg ,(i) 

(2k - r) J ~ n ~ 2 rlB 9 ' {2) 



1,(1) , 1,(2) 

igq 1 '91 



(C.5) 



where B 9 >^ and B 9 >^ 

are B 9 (k,r) evaluated at the poles fe_ and k^\ respectively, 

with 

k 2 k 2 



_(2) 



.(1) = 

2k + 2ki + u 
k 2 -k 2 



2(k + -k 1+ ) 2Jfei+(l-w) 



& k l = 2k+k^ - k = 0, 
h 2 

& k 2 - 



(l-u) 



For the evaluation of B 9 (k,r) we have 
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B«(k,r) = 7m( h~ ftlvD^k - r)V$ a D a ,(k) 

= 7/x(A- /k)luV^{k,k - r) = 7m (/; 1 - jk)^ v {g vtx n-kn-(k - r) 

—n-(k — r)k u n^ — n-n u (k — r) M + — r)n v n^) 
= n-kn-ik-r)^^/;!- jk^^ - n-k{jk- f){jk x - A 

— n-(k — r) jh{jki— Jk) fc + k-(k — r) jh{jk\ — Jk) fx 
= —2n-kn-(k — r)(^i— jk) — n-k{jk jk\— / ^1+ jk — k 2 ) fi 

—n-(k — r) jh{Jki jk — k 2 ) + k-(k — r) jh{jk\— jk) jh 
= -2k + (k - r)+(jki- jk) - k + (/c jk x - / fa+ / jk - k 2 )^ + 

-{k - r) +7+ ( A jk - k 2 ) + (k 2 - fc-r) 7+ (A- $7+ • (C6) 

As before we are interested in picking up the leading logarithmic divergence terms 
which are again quadratic and come from terms proportional to k_ (and k 2 for the pole 
k^). Since 7 OT is related to the gluon-quark transition in evolution, the contributing 
terms will also have to have the structure of a quark vertex, that is, be proportional 
to 7 + as well. We obtain 

B°(k,r) = -2k + (k-r) + (-k_ 1+ )-k + ( 1+1 _ 1+ (k_k 1+ )-0-0-k 2 1+ ) 

-(k - r) + ( 7+ 7_7 + A; 1+ A;_ - k 2 ^ + ) + (k 2 - fc_r+)7+7_7 + (A;i - k)+ 
= 2k + (k — r)+A;_7 + — 2^ + k+k_ki + + r y + k + k 2 — 2(k — r) + ^ + ki + k_ 
+(k - r) + k 2 1+ + 2k 2 1+ (k 1 - k)+ - 2(h - k) + k_r +1+ . (C.7) 



Thus, for the k_ pole we have 
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B 9,(l) = , v pfc + (fc-r) + 2k + k 1+ | Q 2(k-r) + k u 



2k+ 2k+ 2k A 



2(h-k) + r ± \ e 



2k + 

1+ k%k - r) + - fc 1+ - (* - r) + ^± - r+ ih^3± 

127/ U ~ V 1 ~~ U \ 
7 + K fcl+(Co> — V — 1 f ) 

UJ UJ 

^ 1+ k 1+ (uj 2 -2uj) , (Ci 
a; 



while for the fci pole we obtain 

B g,(2) ( 2fc + (fc-r) + 2k+k 1+ _ k+ 2{k-r) + k 1+ 

1+ \ 2(1 -uj)k 1+ 2(1 -uj)k 1+ 1-LU 2(1 -uj)k 1+ 
(k - r)+ 2(A; 1 - fc)+ + 2(£; 1 - £Q + r + \ fc2 



1 — u) 1 — uj 2(1— u)ki + 

l + k 2 



1-LU 



{-uj{k-r) + + k+- k + + (k-r) + - (k-r)+- 2k 1+ (l-u) + r+(l-u)) 



l+k 2 h + / 2 o , o \ 

(a; — 2uj + 2 — v) 

1 — u 

.JL kl+1+ (l + (l-u) 2 -v) . (C.9) 

1 — UJ 



Subsitute these back into ( |C.5|) we have 
-L I ct s Cp 



7, 



.</'/ 



2tt 

+ 

Ol s Cf 

2tt 



4r 
l 2- 



" du;u; cu 2 — 2cj 



1/ 



da,(l - uj) 1 + (1 - uj) 2 - v_ ki+k j 2{2u _ u)M " 



1 - v 



(CIO) 
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upon which we have recovered the result listed in ( 4.53 ), that is, after factoring out 
the logarithmic divergence in the transverse momentum, we obtain 



^(+) 
igq 



2~ix 



Cpj+kt 1 ^ [ dcu(2u-v) M ( Q(v-uo)-{uj 2 -2uu) 
Jo \ v 

-0(^)-^(l + (l-c) 2 -z,)) . 

l-v J 



(c.n; 
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Figures 
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Figure 1: Diffractive Vector Meson Production 
This is the diagrammatic representation of diffractive vector meson production (Part 
(a)). Part (b) is the two-gluon exchange diagram that has the dominate contribution 
when the double leading logarithmic approximation is taken. 
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Figure 2: Deeply Inelastic Scattering 
This is the diagrammatic representation of the Deeply Inelastic Scattering process. 
This is a neutral current process via one photon exchange. The initial proton is \p, a), 
the initial and final electron are k, A and k', A', respectively, and n labels collectively 
the different final states of the proton. 
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Figure 3: Forward Amplitude and Structure Function 
This is the graphic representation of the DIS forward amplitude (a) and the DIS 
structure function (b). The difference lies in the cut (the dotted line) in the middle 
of (b) which puts all the intermediate particles it passes through on the mass shell. 
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Figure 4: Analyticity of the Forward Amplitude 
This figure shows the analytic properties of the forward amplitude on the complex uj 
plane. The crosses are poles and the bold regions on the real axis are the branching 
cuts. The apmlitude has a convergence circle of unit radius, c is the original contour 
used in the dispersion integral. It lies completely within the unit circle and thus 
encloses no singularity. 
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Figure 5: Dispersion Integral Contour 
This figure shows the contour d of the dispersion integral of the forward amplitude. 
It is a continuous distortion from the c contour of figure d does not cross any cut 
nor does it enclose any singularity on the complex uj plane. 
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(1,1) 



x 



Figure 6: Integration Region and Change of Variables 



This figure illustrates the change of variables and integration limits in section p.. 4.2 



where we show the equivalence between the DGLAP equation and the RGE from 
OPE analysis. The shaded area is the integration region. 
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Figure 7: Forward Operator Product Expansion 
The operator product of the forward amplitude to the lowest order. The upper part 
of the righ hand side is the Wilson coefficient. The lower part is the reduced matrix 
element. We have only shown the leading quark-quark transition. 
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Figure 8: The Lowest Order Quark Distribution Function 
This is the quark distribution function pf(x) where we have put in explicitly the 
lowest order quark vertex in the light cone gauge. The dotted line is the cut that 
puts all the intermediate lines it passes through on the mass shell. 
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Figure 9: The Lowest Order Quark-Quark Diagram 
This is the quark-quark diagram of the first order radiative correction to the quark 
distribution function, a, b are color indices and i labels the on-shell gluon. 
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Figure 10: Quark Self Energy Diagram 
These are the two quark self energy diagrams that are needed to be added to figure |9] 
to cancel infrared divergence at the end points of integration. 
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Figure 11: The Quark-Gluon Mixing Graph 
This is the graph that generates the anomalous dimension for quark to gluon mixing 
in forward DIS. 
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Figure 12: Non-forward Scattering Amplitude 
This is the unequal mass double virtual Compton scattering diagram. Both q and q' 
are virtual photons and they are of different invariant masses. There is a non-zero 
four-momentum transfer r from the incoming virtual photon to the target proton. 
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Figure 13: Symmetry of the Non-forward Scattering Amplitude 
Because the two diagrams shown are indeed identical after the momentum labeling 
change /i <-> v, q <-> —q', the non-forward amplitude is invariant under the same 
transformation. 
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Figure 14: Handbag Diagrams in QCD-Parton Picture 
These are the so-called "handbag" diagrams contributing to the non-forward am- 
plitude in a QCD-Parton picture. We have explicitly included the cross diagram. 
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Figure 15: Triangle Graphs for the Anomalous Dimensions 
These are the four lowest order "triangle" diagrams contributing to the anomalous 
dimensions. The momentum labels of the four are the same (only the quark-quark 
diagram shows them explicitly). 



160 




Figure 16: The Quark-Quark Diagram of the Anomalous Dimensions 
This is the detailed version of the first triangle graph, the quark-quark diagram, a 
and P are Lorentz indices; a, a' , b and b' are color indices for the quarks while % labels 
the qluon line. The difference between this diagram and Figure || in the forward case, 
besides the non-forwardness, is that we are not putting the gluon on-shell. 
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Figure 17: The Poles and Contours of k_ Integration 
The figures show, on the k- complex plane, the poles and different contours of fc_ 
integration in the four regions of its value (see section [4.2. 1| ). The dots are the pole 
positions with their distances to the real axis exaggerated. 
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Figure 18: The Tree Level Born Diagrams of the Non-forward Amplitude 
These are the tree level (Born) terms contributing to the non-forward amplitude 
from which we can calculate the lowest order Wilson coefficients. We have explicitly 
included the cross term. 



163 




Figure 19: The Gluon-Quark Diagram of the Anomalous Dimensions 
This is the detailed version of the second triangle graph, the gluon-quark diagram. 
a, P, \i and v are Lorentz indices; a, a', and b are color indices for the quarks while i 
labels the gluon lines. 
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Figure 20: The Quark-Gluon Diagram of the Anomalous Dimensions 
This is the detailed version of the third triangle graph, the quark-gluon diagram, a, 
(3 are Lorentz indices; a, b are color indices for the quarks while i and i' label the 
gluon lines. 
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Figure 21: The Gluon-Gluon Diagram of the Anomalous Dimensions 
This is the detailed version of the last triangle graph-the gluon-gluon diagram, a, 
a', j3, j3', 7, 5, ji and v are Lorentz indices; i, j, k and k' label the gluon lines. 
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Figure 22: The Gluon Self- Energy Diagram 
These are the gluon self-energy diagrams that are needed to cancel the end-point 
divergences in the LC calculation. We need to include both the gluon loop and the 
quark loop contributions. 



